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ABSTRACT OF THE DISSERTATION 


Physics of Nonmagnetic Relativistic Thermal Plasmas 

by 

Charles Dennison Dermer 
Doctor of Philosophy in Physics 
University of California, San Diego, 1984 
Professor Robert J. Gould, Chairman 

A detailed treatment of the kinematics of relativistic systems 
of particles and photons is presented. In the case of a relativistic 
Maxwell- Boltzmann distribution of particles, the reaction rate and 
luminosity are written as single integrals over the invariant cross 
section, and the production spectrum is written as a double integral 
over the cross section differential in the energy of the produced 
particles (or photons) in the center -of -momentum system of two 
colliding particles. The result is valid for all temperatures and for 
the general case when the plasma contains different mass particles. 

The results are applied to the calculation of the annihilation 
spectrum of a thermal electron-positron plasma, confirming pre- 
vious numerical and analytic results. Relativistic thermal electron- 
ion and electron- electron bremsstrahlung are calculated exactly to 
lowest order, and relativistic thermal electron-positron 


ix 



bremsstrahlung is calculated in an approximate way, by interpolating 
between the modified Bethe-Heitler formula and the exact lowest 
order electron-electron bremsstrahlung spectrum. An approximate 
treatment of relativistic Comptonization is developed. 

The question of thermalization of a relativistic plasma is 
considered. A formula for the energy loss or exchange rate from 
the interaction of two relativistic Maxwell- Boltzmann plasmas at 
different temperatures is derived. It is concluded that at tempera- 
tures greater than about 1 MeV, full thermalization of an electron- 
proton plasma is not possible because of the dominance of brems- 
strahlung losses. 

Application to a stable, uniform, nonmagnetic relativistic 
thermal plasma is made. Approximate expressions for dominant 
rates, taking account of the effects of Comptonization, are used to 
model the system at extreme relativistic temperatures and at sub- 
relativistic temperatures in the limit of large optical depth. 
Equilibrium pair densities, critical temperatures, and emergent 
spectra are calculated. Comparison is made with other studies. 
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I. INTRODUCTION 


The existence of relativistic plasmas is supposed to account 
for many astrophysical observations made during the past two 
decades. These observations include gamma ray bursts, jet struc- 
tures observed near some extragalactic radio sources or stars 
(SS433), and the central energy sources of active galactic nuclei and 
quasi-stellar objects including, perhaps, the center of our own 
galaxy. One might also include indirect observations relating to that 
epoch of the early universe when nucleosynthesis and neutrino de- 
coupling were occurring, according to standard cosmology. In order 
to characterize the physical conditions relating to such observations, 
a systematic study of relativistic gases is required. 

A relativistic plasma, or gas, is a system of particles and 
photons in which the average particle and photon energy is suffi- 
ciently great that particle -creating reactions occur and influence 
the observed properties of the gas. The created particles are pro- 
duced at the expense of excess kinetic energy in collisions between 
two particles in the plasma, or as a result of particle production 
through photon- particle or photon-photon interactions. Since the 
produced particles may in turn participate in particle -creating 
reactions, the description of the system requires a self-consistent 
treatment which is, in general, nonlinear. 

The aim of this work is to examine, by means of theoretical 
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analysis of relativistic particle gases, processes in relativistic 
plasmas, in particular, investigation is made of processes occur- 
ring in the "tr an s relativistic" regime where the average particle 
kinetic energy is of the same magnitude as the particle rest mass 
energy. Because cross sections for particle or photon production 
can be extremely complicated in this regime, asymptotic forms may 
not be used and a simplified, general treatment is essential. Of 
course, the results must properly approach limiting forms in the 
nonrelativistic and extreme relativistic limits. 

Conditions of thermal equilibrium cannot be assumed in the 
study of relativistic gases. Indeed, the predominance of certain 
energy-loss processes can prevent the formation of a thermal dis- 
tribution of the particle momenta (Gould 1982a; Stepney 1983), as we 
shall examine in more detail later. The term relativistic gas will 
be used interchangeably with the term relativistic plasma; though the 
system will be in a state consisting of independent charged particles, 
the possibility of collective modes of excitation will not be treated in 
this study, as we consider a nonmagnetic system. 

The examination of the properties of a relativistic gas begins 
with an elaboration of the various processes that contribute to the 
behavior of the gas. The processes which are most important de- 
pend on the nature of the problem under consideration. It is then 
essential to calculate energy loss rates and time scales in order to 
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determine whether the particles of the gas can be thermalized, as 
the form of the particle distribution function defines particle produc- 
tion and annihilation rates and photon spectral emissivities. 

It will then prove useful in a self-consistent analysis of the 
nonmagnetic relativistic plasma to consider ranges of system 
parameters such as optical depth, proton number density, and tem- 
perature, where expressions for rates assume convenient analytic 
forms. In addition, an exact treatment of the Comptonization of 
internally generated photons by a relativistic particle plasma is a 
formidable task, so approximate techniques are necessary. These 
results can be compared to a gross numerical simulation. 

The plan of this dissertation is as follows: A theoretical 
formalism is developed in Chapter II for the calculation of rate 
processes and timescales in relativistic gases. A general formula 
for the production spectrum from an isotropic relativistic Maxwell - 
Boltzmann gas is derived, and a treatment of the Compton power and 
Compton scattered spectrum is developed. Timescales for reaction 
and energy loss processes are also examined. Application of the 
formalism is made in Chapter III, including the calculation of annihi- 
lation spectra and bremsstrahlung from a relativistic Maxwell - 
Boltzmann plasma, calculations of absorption and opacity, and a 
compilation of invariant cross sections useful for the analysis of 
problems in high energy astrophysics. 
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Physical applications of the preceding results are developed 
in Chapter IV. A uniform nonmagnetic relativistic plasma is ana- 
lyzed in various temperature and optical depth regimes, asymptotic 
expressions for the dominant rates are derived, and an approximate 
treatment of Comptonization is developed. Discussion of the prin- 
ciple features of this system follows , including comparison with 
other work. Results and conclusions are set forth in Chapter V. 



n. THEORETICAL TREATMENT OF PROCESSES IN 
A RELATIVISTIC PLASMA 


A. Introduction 

Historically, the theory of the relativistic gas begins with the 
derivation of the relativistic generalizations of equilibrium particle 
distribution functions (Juttner 1911a, b; 1928). Later developments in- 
clude generalizations of kinetic equations for relativistic gases. 

Early results are summarized in a monograph by Synge (1957), and 
details of the kinetic theory of relativistic gases may be found in the 
book by de Groot, van Leeuwen, and van Weert (1980). 

In this chapter, expressions for rate processes and the spec- 
tral emissivity are derived for isotropic relativistic particle gases. 
The results are specialized to the case of a relativistic Maxwell- 
Boltzmann (MB) distribution. Considerable simplification of the 
production spectrum of a MB gas is achieved in this case. Expres- 
sions for reaction and energy-loss timescales are also presented. 
Finally, a formalism is developed for calculating the Compton power 
and approximate Compton scattered spectra for a gas containing 
photons and electrons. 

B. Relativistic Kinematics 

Consider a relativistic gas of particles. The relativistic 
reaction rate r is defined as the number of collisions per unit time 
per unit volume and is an invariant quantity. The reaction rate for 


5 



6 


two interacting distributions of particles is given by 


r 


( 1 + 6 i 2 ) 


-1 





(1) 


(Landau and Lifshitz 1962; Gould 1971). In this expression, com- 
pletely general, the integration extends over all momenta of the 
interacting distributions of particles. The number of particles per 
unit interval of momentum in the observer's system of reference 
(LS), or the spectral density, is denoted by dm for particles of the 
type i. The factor (1 +6 _) * corrects for double counting if the 

J. u 

particles are identical (self-interacting) and da is the invariant 
cross section. 

The speed of light c is set equal to 1 in equation (1) and the 
momentum p. , energy E., and mass m. are written in energy 
units throughout. In this case (3. = p./E. and {3 , the invariant 

AM 1 1 1 r 

relative velocity, is defined as 




, ,2 2 2 

( tt. • tt_) - m. m_ 

l 2 12 


(T V V 


1/2 


(2) 


where the four-momentum of particle i is given by tt. = (E . , p.) 

and the metric is (+,-,-,-). The variable 8 represents the 

r 

absolute value of the velocity of one particle in the frame of refer- 


ence in which the other particle is at rest. The factor (1-8 • |3 ) 

■MA 1 

accounts for relativistic changes in time and particle density in 
transforming from one reference frame to another. 
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In the general case of two distinct nonisotropic interacting 
gases, equation (1) is a 6-fold integral and admits no simplification. 
Assuming, however, that the interacting gases are isotropic, equa- 
tion (1) may be written as an integration over three quantities neces- 
sary to uniquely characterize a collision between two different 


particles. These are conventionally taken to be the momenta p^ 
and p , and 9, the angle between the direction of the two momenta 

Lt 

in the LS. Thus, cos 0 = p • p^/p.p-,. In terms of these variables, 


2 03 

8tt A A_ r 


>2 P 2 £ 2 (P 2 } d(GOS 9) 


X [ct(p 1 ,P 2 ,cos 9) . (1 - ^ cos 9) • P r (P 1 »P 2 « cos 9)] 


The fxmctions f.(p.) describe the momentum distributions of the 
l i 

interacting particles and are related to dn . and the normalization 
coefficient A^ through the expression 






( 4 ) 


n. is the volume density of particles of the type i . 

Ordinarily, equation (3) is used to calculate reaction rates 
as it stands. It is reasonable to suppose that equation (3) can be 
simplified by a judicious choice of variables. In particular, one 
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such variable should be the argument of the invariant cross section. 
In this case, the integration over this variable can be performed 
last, and the generally complicated dependence of the cross section 
can be avoided as long as possible. 

To this end, it is convenient to choose the invariant Lorentz 
factor associated with the relative velocity j3^ , defined in equa- 

tion (2), as this variable. In terms of p , p_ , and Q, 

6 1 ^ 


2 - 1/2 V "2 E 1 E 2 * P 1 P 2 COS 9 

v r = (i - e r > 


m i m 2 


m i m 2 


(5) 


For the other two variables, introduce v and the angle uu, where 
Y c is defined as the Lorentz factor corresponding to the velocity of 
the center -of -momentum of the two interacting particles with re- 
spect to the LS. It can be written as 


Y 


c 



+ E, 


')] 


1/2 


(6) 


in terms of the invariant quantity 

2 2 2 

s = (rr + tt 2 ) =m 1 + m 2 + 2m ] m z y r , (7) 

which describes the strength of the collision between the two par- 
ticles. The variable v is a natural extension of the c enter -of - 

' c 

momentum variable used in nonrelativistic studies of rate 


processes. 
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The angle tu is defined through the relation 


sm a) = 


P 1 P 2 sin 0 


mm 8 y 8 Y 
1 2 r T r c < 


(8) 


The significance of uj is as follows: If we transform to the center- 

of -momentum system (CM) of the two colliding particles and call the 

transformed momentum of particle 1 p* , then uu represents the 

angle p ' makes with respect to the direction of transformation 

from the LS to the CM (see Figures la and lb). ' 

The Jacobian J for the transformation from the variables 

p , p , cos 8 to y , y » cos <i> is found to be 
12 1 r 1 c 


'P x » P 2 . cos 9 


Y r » Y c » cos 0), 


2 2 

E E m m g y I 3 Y 
1 2 1 2 K r T r r c T < 

2 2 

Pi P 2 


(9) 


Changing to this new set of variables in equation (3), it is convenient 
to rewrite equation (3) in a general fashion so as to incorporate ex- 
pressions for different rate processes. Letting P(y) stand for a 
rate process of the relativistic gas for the parameter value(s) y 
associated with the weighting. p.(y; Y r ' Y c » ^)» and defining 
u = cos uu , then 


2 3 3 ® 

8tt A A rn m 


P(y) = 


l A 2 m i m 2 f 2 f • 2 
TT^ J d Y r W r - 1) J d Y c (Y c - 1) 


1/2 


( 10 ) 


/ du f x [p 1 ( Yr , Y c » u )3f 2 [p 2 (y r » Y C » U )I P (y; Y r » Y c ' u ) • 

-1 


X 
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Figure 1. Kinematic quantities associated with relativistic 

transformations. Transformation from the LS to the CM is 

effected by a Lorentz transformation by momentum p de- 

*+* c 

fined in terms of the momenta Ip I and ]p | of the inter- 
acting particles and the angle 9 between them. Transforma- 
tion to the C- system is effected by an additional counter- 
clockwise rotation of the axes by an angle uu. The terms 

❖ 

p, p 7 and p represent the momentum of a produced par- 
tide as seen in the three frames of reference. The coordi- 
nate system employed in the three frames (illustrated only 
in the C-system) consists of specifying the absolute value of 
p, the cosine of the angle ip that p makes with the z-axis, 
and the angle 0 between the x-axis and the projection of the 
momentum p on the x-y plane. 
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For example, P(y) = r, the reaction rate, when P (y ; y y , u) = 

r c 

CT (y^). The luminosity L is the total energy produced in particles 

or photons per unit volume per unit time as measured in the LS. If 

the energy liberated per interaction as recorded in the LS is denoted 

by e(Y , Y , u), then P(y) = £ when p(y; Y , Y , u) = 
r c re 

cr (Y r ) e (Y r , Y c > u). In like manner, the production spectrum S(E) 

(number of photons produced per unit volume per unit time per unit 

energy interval dE) is P(y) = S(E) when p (y; Y , Y , u) = 

" ' r c 

do(E; y y , u)/dE. This last quantity is the cross section differ- 
r c 

ential in the LS energy of the produced particles or photons in an 
interaction characterized by the reaction parameters y Y c » 
and u. 


These rate processes are closely related to one another, and 
consistency of results requires that the following relationships hold 
for the reaction rate r, the luminosity £, and the production spec- 
trum S(E): 


and 


r = / dE S(E) (11) 

m. 

J 



dE E S (E ) , 


( 12 ) 


where m. is the rest mass of the particle produced in collisions 
between particles in the gas (m. = 0 if photon production is 
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considered). 

We now specialize to the case of collisions in a particle 
plasma described by a relativistic MB distribution in order to de- 
rive an expression for the reaction rate. The MB distribution for 
relativistic particles of the type i is 

f. (p.) = exp [-b(pf + mf )* ] , (13) 

11 l i 


defining b = 1/k^T and noting from relation (4) that 

n - M • 

A 11 

i 3 

4 TTmf KfP .) 

l 2 i 


(14) 


where p . = bm. and K (x) is the modified Bessel function of 
1 1 v 

order v (Synge 1957). Assuming that both interacting gases have 
MB form at the same temperature, equation (6) implies that 

fl(Pl) f 2 (p 2 ) = exp t-b-Y c [s(Y r )] 1/,Z } . (15) 


Substituting equation (15) into equation (10) for the case of the reac- 
tion rate gives 


r = 


n, n 2 u ; u 2 




a 

/ 


d y. 


(Y r Z - DKjtblslY^] 172 } 

tb[»(Y r )] 1/2 } 


a (y ) 

r 


( 16 ) 


agreeing with the result of Weaver (1976). 
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C. Secondary Particle Production Kinematics 

The preceding considerations show how the reaction rate and 
luminosity are related to the spectral emissivity. In order to derive 
the spectral emissivity, it is necessary to calculate the LS cross 
section differential in energy of the produced particles or photons 
for an arbitrary collision between two particles. In terms of the 
variables y r » Y , and u, it is clear that all particle reactions 
characterized by the value y^ are kinematically equivalent in the 
C-system, the reference frame rotated so that the directions of the 
incident particles are along the z'“'-axis of the C-system (see Fig- 
ure lc). It is then only a matter of performing the appropriate 
Lorentz transformation from the C-system to the LS, defined in 
terms of u = cos (l) and y^, to get the cross section differential 
in the LS energy of the produced particles or photons. 

The appropriate transformation consists of two successive 
Lorentz transformations: a clockwise rotation of the coordinate 
axes of the C-system by an angle u) followed by a boost of the CM 
by velocity - (3 along the x' axis. The equations of transforma- 
tion are easily derived, and it will be useful for later work to dis- 
play the relevant coordinate transformation equations. Quantities in 
the C-system are labelled by asterisks. Define a coordinate system 
in the LS and C-system by the variables (p, cos 4 1 , 0) and 

«T* *T* 

(p"', cos vjj , 0 ), respectively, as illustrated in Figure 1. A 
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produced particle whose energy is E , momentum is p^ , and 

«JU 

'I' 

production angles are cos ^ and 0 has energy 


E = y [E + P^(p cos 4 cos uu - p sin 4 cos 0 sin to ) ] (17) 


in the LS. The reverse transformations are given by 


E = y (E-P pcosiH, and 
c c 


(18a) 


p cos 4 = y (p cos 4* - 3 E)cos to + p sin 4* cos 0 sin u). (18b) 


To determine the cross section differential in energy in the 
LS, it is best to use invariant production cross sections, as the 
question of the Jacobian of the transformation is thereby avoided. 
The Lorentz invariant cross section for the production of particle j 
in the reaction 1 + 2 j + X is given by 


d 3 a 


dp 


* 

a a . 
L 


*T> -1' 'I 

p dE d(cos 4 ) d 0 


* 


F (E 
J 


*3 


Y ) 

r 


(19) 


X represents all final state particles other than those of the type j. 

* l i' 

The invariant cross section F. represents the probability that the 
produced particle j will be deposited in the phase space volume 

3 ❖ 

d p /E in a collision characterized by the reaction strength s, or 
equivalently, Y r » since the two are linearly related through equa- 
tion (7) (Berger 1975). Asterisks in the definition of the Lorentz 

*T* 

invariant cross section in equation (19) serve to remind that F. is 
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most conveniently written in C- system quantities, and in fact experi- 
mental results are usually presented in just this coordinate system, 
though the LS and CM representations are completely equivalent. 

The invariant cross section is normalized through the relation 


(Y r » 


TOT' y r 


(y ) = 


/ 


3 Jj 

d 3 P. 


F. (E 


Y ) 

r 


( 20 ) 


with the integral extending over all momenta. In this expression 

( avera g e multiplicity of the produced particle j and 

CT^^^Cy ) is the total cross section of the reaction 1+2 character- 
TOT r 

ized by the collision strength Y r » though other conventions are pos- 
sible (Likhoded and Shlyapnikov 1978). 

The energy spectrum in any frame related by a coordinate 
transformation, whether it be a rotation or boost or some combina- 
tion of the two, is given by 


da(E; y 


u) 


dE 


= (E 


2 , 1/2 

m. ) 

J 


/ d0 / 


d(cos i\>) F. (E , 4» ;y ) (21) 

J r 


-i' 

with E and ^ written in terms of E, 4* , and 0 according to the 
transformation equations (18). Azimuthal symmetry in the produc- 
tion of j has been assumed in equation (21), and the limits of inte- 
gration over cos 4 1 and 0 remain to be specified. 

■A. 

Define ? = as the maximum possible energy that 

the produced particle j can carry off in a reaction characterized by 
the value Y p • Relativistic mechanics can be used to explicitly 
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2 2 1/2 

evaluate this quantity, giving 5 = (s - m^ + rm )/2s , where m^. 

is the smallest possible value of the mass of all remaining particles 
so that the reaction 1 + 2 -» j +X may proceed. The limits of inte- 
gration are then implied by the transformation equations (18). For 
fixed energy E, cos ^ extends from cos = (y - | )/|3 ^ y P to 
cos 4* = 1. The integration is over all 0 as equation (18a) is inde- 
pendent of 0 (Taylor et al. 1976; Tan and Ng 1981). 

Employing the preceding considerations with equations (19) 
and (21), the following identity may be written: 


da (E; y y , u) 
r c 

dE 


= P 


2tt 1 ^ 

/ d0 / d(cosi(j) / dE Q 

0 cos *'m. 

M j 


2 *.U 

*T* *T» 'I' 

d a. (E cos 4> ; y ) 
J 0 _£_ 

dE^ d(cos 4 1 ) 


6(E -E Q ) 
2TT Pn 


( 22 ) 


«t* 'i' 

with E and cos i ]t given by equation (18). For the problem of the 
production spectrum S(E) of a relativistic isotropic MB gas, equa- 
tions (10), (15), and (22) give 


S(E) =- 


n lVl U 2 


2 (l+6 12 )K 2 ( Ml )K 2 ( M 


J <V V r ‘ l) JdY^Y 2 -!) 1 


/2 


1 /2 

exp { -bY c t s (Y r )] 1 
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x 


P 

2tt 


§ * . 1 2tt 1 

[ —T / du / [ d(cc 

m. P 0 - 1 0 cos 4L , 

1 M 


d(cos 4 1 ) 6( E - E ) 


a. o- 


d a . (E cos <\> ; y ) 
3 0 £ 

*A» ^ 

dE^ d(cos 4< ) 


(23) 


Call the term in brackets I. Performing the integral over cos 4* 
using equation (18a) in the delta function and equation (18b) for 
cos ^ , one finds that 


I = 


A * 1 2tt 

2TTB 1 y / — p f du f d« 

C C y (E - (3 p) P 0 - 1 ^0 

c c 


? ❖ * * 2 1/2 2 1/2 
d a. [E Q) cos 4 = pu + (1 - p ) (1-u) cos 0 ; y 

T «T* 

dE Q d cos 4 j 


(24) 


where p is defined as (E - y E )/B y p„ and is bound by the 

c U ecu 

limits -1 and 1 due to the previously established limits on cos 4*. 

The last two integrals in equation (24) represent an integra- 
tion over all solid angles in u - 0 space. Identifying cos 9^ with 
p, we may substitute cos y for cos 9^ cos uo + sin 9^ sin U) cos 0 
using a standard trigonometric identity. The limits of integration 
remain as before -- over all solid angles -- but now over an inte- 
grand possessing azimuthal symmetry. The expression becomes 
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, r dE r 

I= — / ~r / dfcos *> 

c C V P 0 J -1 


d a. (E cos x ; y ) 
j 0 r 

• j , 

dE Q d(cos x) 


(25) 


The last integral in equation (25) is just the production cross section 

«)y 

integrated over the angle X (° r equivalently, ^ ) in the C- system 

A s'- J- 

'I- T* "I" 

and will be denoted as dCJ (E 0 ; y )/dE . (This result reflects the 

Or 0 

uniformity in cos 0 ) of the phase space factor in equation (10) and 
the independence of the product MB distribution, equation (15), of 
cos (JO . ) 


The production spectrum (23) becomes, reversing the order 


of integration over E^ and Y c > 


S(E) 


1 n l n 2 U l P 2 / 2 

2 <1+8 12 >K 2 <Ui>K 2 (h 2 )/ 



s'- 

da (E Q ; y r ) 




, r / n 1 / 2 

_ by c [ s (“Y r ) ] 


(26) 


The limits of integration are implies by a consideration of those 

•A. 


values of y which contribute for fixed values of E and E, For 
c U 

material particles, y = (EE ± pp )/m. , whereas for photon pro- 

C U U J 

2 *2 * + 

duction the limits are y = (E + E„ )/ZEE. and y -» 00 (with, 

c 0 0c 


s'- O- 


of course, = 0 and p^ = E^ in the integration over E^ ). 

Using these limits, the integral over y^ is trivially per- 
formed and equation (26) gives the general formula for the 
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calculation of the production spectrum of particles or photons from 
interactions in a MB gas at all temperatures, from the nonrelativis- 
tic to the extreme relativistic limits. It is written in terms of the 
cross section differential of the produced particles (or photons) in 
the C- system of the collision. As most cross sections are calculated 
or measured in terms of these variables, equation (26) is a con- 
venient formula for analysis. Clearly no further simplification can 
be made without referring to a specific problem. Applications of 
equation (26) will be deferred to the next chapter. 

From equation (26), the reaction rate and luminosity can be 
derived through relations (11) and (12). It is straightforward to 
integrate the production spectrum (26) over energy to recover the 
reaction rate (16). To obtain a corresponding expression for the 
luminosity, define 


e 


(y ) 

r 


(y ) 

r 



da (E 


0 


y ) 

r 



(27) 


as the average energy, weighted by the cross section, of particles or 
photons released in the CM of a collision characterized by the reac- 
tion strength y^ . Then the luminosity is 


t = 


n n M- M> 
12 12 


(l+6 12 )K 2 (Hi) K 2 (H 2 ) 


0 

/ 


dy_ 


(Y^-l)K 2 tb[.(Y r )] 1/2 ) 

fb[s(Y r )] 1/2 ) 


e (y )a (y ) 
r r 


( 28 ) 
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D, Annihilation and Energy Loss Timescales 

The derivations of equation (26), the spectral emissivity, as 
well as the expressions for the reaction rate and luminosity, equa- 
tions (16) and (28), are performed under the assumption that the 
particles are described by an isotropic energy distribution of MB 
form. This assumption considerably simplifies the results. Wheth- 
er, or under what conditions, this assumption is warranted will be 
the concern of this section. 

It should be remarked at this point that under many condi- 
tions of astrophysical interest, the state of complete thermodynamic 
equilibrium, in which the photons are of Bose-Einstein form and the 
electrons and positrons possess a Fermi- Dirac distribution at zero 
chemical potential, cannot be achieved. Yet the particles can be 
thermalized, that is obtain MB form (which is just the Fermi- Dirac 
distribution with large negative chemical potential). The establish- 
ment of a thermal distribution of particles requires that elastic 
scattering processes dominate over inelastic processes. The proc- 
esses that drive a gas composed of electrons, positrons, and ions 
to thermal form are the elastic Coulomb scattering processes, 
namely Miller scattering between two electrons or two positrons, 
Bhabha scattering between an electron and a positron, and Ruther- 
ford and Mott scattering between ions and an ion and electron (or 
positron), respectively. Bremsstrahlung or synchrotron emission 
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are typical energy loss mechanisms that can prevent the establish- 
ment of a thermal distribution. 

There is a considerable literature dealing with the thermal- 
ization of particles at subrelativistic temperatures. Basic results 
can be found in Spitzer's book (1962). A systematic study of the 
thermalization of particles in a relativistic gas has only recently 
been considered. Gould (1981b) treats the thermalization rates of a 
highly relativistic particle in the presence of a MB gas at tempera- 
tures greater than 0.5 MeV. Since the high energy particle is in the 
tail of the Maxwellian of the relativistic gas, which is in fact that 
portion of the distribution function requiring longest to relax, the 
timescales for energy loss accurately represent the time of a 
particle gas to relax to thermal form. 

In subsequent papers (Gould 1982a,b, c), comparisons of the 
earlier results are made with the energy loss timescales due to 
various inelastic processes in relativistic electron and nucleonic 
gases, and it is shown under what conditions a thermal distribution 
can and cannot be obtained. In particular, above a certain tempera- 
ture, a relativistic electron gas cannot be thermalized because of the 
predominance of bremsstrahlung losses. Moreover, a sufficiently 
strong magnetic field can prevent thermalization through synchrotron 
energy losses at even lower temperatures. 

The approach to be taken here is somewhat different. First, 



a general expression for the reaction rate of a single test particle in 
the presence of a MB gas will be derived. This result will provide, 
for example, the timescale for annihilation of a positron in an elec- 
tron gas, or the timescale for an energy loss process in which a 
large amount of energy is lost in a single collision as occurs, for 
instance, in pair production through particle-particle collision. 
Second, an exact expression for the energy loss of a test particle in 
the presence of a relativistic MB gas is derived. Then the timescale 
for energy exchange between two MB gases will be determined. This 
result is a generalization of an expression derived by Stepney (1983), 
and is valid for two interacting MB gases at any two temperatures. 

Suppose a particle of momentum p° interacts with the 


particles of a relativistic gas described by the momentum distribu- 
tion f^(pj). The reaction rate is given from equation (10) by 


2 3 3 00 oo 1 

8TT A 1 A 2 m 1 m 2 [ r 2 l/Z f 

r = i + 6,1 / dY r ( V 1} / dY c (Y c" 1) / duf l ( Pl 

, 1 


)f 2 (p 2 )a(Y r ) 


(29) 


Since the reaction rate is the number of collisions per unit volume 

per unit time, the average time for the particle of momentum p° to 

survive between collisions is T = (r/n ) \ where n_ is the 

r u Cd 

o 

density (to be normalized to unity) of the particle of momentum p^ , 
and is closely associated with the absorption coefficient, as will be- 
come clear in Chapter III. 
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The distribution function of particles of the type 1 is now 
taken to be of MB form, given by equation (13), with normalization 
coefficient (14). The momentum distribution function of the particle 
of momentum p° is a delta function. Incorporating the phase space 
factor, it is convenient to take the distribution function as 


f 2 (p 2 ) 


4 TTpJ 


6(p 2 -p 2 ) 


(30) 


From the normalization condition (4), - n^ = 1. 

The reaction rate (29) for an individual test particle becomes 


2 00 00 

n b m m 1 f 2 f 2 

= ' 1 l ■ ■ * — / dy (y -l)cr( Y ) / dy Cy - D 

2K (b m ) o o / r r r / c c 

2 11 A 


1/2 


2 2 1 


L 


-b E 

due 6(E 2" E 2 ) » 


(3D 


where b = 1/k^T refers to the inverse temperature of particle 
distribution 1 and the delta function has been converted to energy 
variables. From a consideration of equations (5), (6), and (7), it is 
possible to express E^ and E^ in terms of the variables y^ , y^ , 
and u. The relations are 


E i = 


m (m y +m.)y 
1 2 r 1 c 


1/2 


mm 8 y P y 
1 2 r r c c 

1/2 


u 


and 


(32a) 
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E 2 


m. (m y + m ) y m m (3 y P Y 
2 1 r 2c 1 2 r r c c 


1/2 


1/2 


u 


(32b) 


Performing the delta function over u and determining what the 
conditions -1 £ u ^ 1 implies for the limits on y we find that, 
after integrating over y , the reaction rate r^ is 


n. 


1 1 
r t ' 2K (Uj) 


P 2 Y 2 1 


d Y P_ y„ 0 (y ) 

r r r r 


-UlY 2 Y r 


W I fJ 2 Y 2 fi r Y r -^ fS 2 Y 2 f 'r Y i 


- e 


(33) 


where = ^.m., , and m 2 y 2 B 2 = p° . 

Equation (33) gives the formula for the annihilation or reac- 
tion rate of a single test particle of momentum p° and mass m 2 in 
the presence of an isotropic MB gas of particles of mass m^ at 
temperature 1/b^. To be consistent with earlier results, it should 
reduce to equation (16) when integrated over a MB distribution of 
particles of type 2, also at temperature 1/b^. To determine the 
reaction rate of two interacting MB gases at different temperatures, 
note that 


' b 2 E 2 


d P 2 (4 TTp 2 2 ) e 2 2 f 2 (p 2 > 


(34) 


where f 2 (p 2 ) is given by equation (30) and the normalization coeffi- 


cient is now given by expression (14). Multiplying equation (33) by 
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(4rrp°^) exp(-b E°) and integrating over dp° according to equa- 

Ci Cd Ld Cd 

tion (34), the reaction rate is 


2(1 + 6 


12 


x dv 2 e 


/ 

n l n 2 M 2 f 

J d \ 0 r \ ° < Y , 



-(^ + w;v r )Y 2 


^l' 0 r Y r P 2 Y 2 ‘ M l' 0 r Y r ®2 V 2 


- e 


(35) 


The remainder of the derivation is similar to the treatment of 


W 


1/2 


eaver (1976). Defining ch9 1 = (p ' + p ' y )/z and ch0 o - y , 

1 2 1 r 2 2 

2 2 

with z = W- ' + M- * + 2P_' p' v , gives Weaver's result for the reac- 
1 2 1 2 'r ’ 

tion rate of two interacting MB gases at different temperatures, 
namely 


n iVj_!±i_ f S 
: 2 <w P K 2 <U 2 ) / Yr 


r = 


< 1+ V K 


(v 2 - Dc(y )K,( Z 1/2 ) 

r 1 


1/2 


(36) 


If the gases are at the same temperatures, p^ = p^ , P^ = U 2 > 
z = bs , and equation (16) is recovered. 

To determine the thermalization and energy loss timescales 
for various processes in a relativistic gas, consider the expression 
for the change in the LS energy of a particle before and after a colli- 
sion. The average energy change can be defined as 
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<AE> = 


a( Y ) 

r 


,3 ; 

d p 


dV 


dp 


*3 


AE 


(37) 


3 -\'2 * ❖ Y 

where dp = p • dp d cos i(j d$ , employing the coordinate system 

a*/ 

illustrated in Figure 1. The cross section for scattering into the 

3 * 3 ❖ 3 * 

phase volume dp is given by (da /dp) and A E is defined as 

/W AV 

the change in the LS particle energy as the result of being scattered 

3 * 

into the phase volume dp in the C- system. From equation (17), 

<w 

the energy change is just 


AE - y (E - E^ ) + 0^ y c [(p cos 4* - p. ) cos U) - p sin 4* cos 6 sin uu], 

(38) 

where p. and E. refer to the particle momentum and energy prior 

A 

to scattering. The value of p. is in fact given by 


p (y ) = (m m R y )/[s(y )] 

i r i r r r 


1/2 


(39) 


The momentum and energy of the particle after scattering are de- 

* * 

noted by p and E , respectively. 

The timescale for energy loss by a system, whether it be a 
gas or a test particle, is defined as T *= £ /e , where e is the 

i. K K 

kinetic energy density of the system. If the system loses energy by 
inelastic scattering through the production of particles or photons, 
then i has been determined for the case of two MB gases at the 
same temperature, namely equation (28). If the scattering is 
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elastic, or if the energy exchange rate of two MB gases at different 
temperatures is being considered, it will be necessary to proceed 
differently. 

Consider first a test particle that loses LS energy through 
collisions with particles of an isotropic MB gas. If the test particle 
has momentum p^ , its momentum distribution is described by equa- 
tion (30), and the rate of change of the LS energy of the test particle 
is given through equations (10) and (38), after integrating over the 
delta function, as 


b E° 

l = j 1 2 fay e y B 1/2 f dE* f d cos W f) 

\ r F r J, U dcos J 


(40) 


x 


/: 


. 1/2 
-b 8 v 


dV V e 
c c 


. * * 

(E -E. _)+p 
l, 2 c 


* * 

(p CO 8 ^ + p 


l) 


The subscript "t" again refers to the case of a test particle, and the 

maximum and minimum values of the cosine of the scattering angle 

in the C-system, cos ^ , are denoted by 1 and 2, respectively, and 

Y* = (E° E* ±p° p! )/ml . The quantity E* = (p? 2 + mV /2 
c 2 i, 2 2 i, 2 2 i,2 i 2 

and is the energy of the test particle in the CM, noting, from equa- 

tions (32b) and (39), that E = y (E. - 6 p. u). 

2 c i, 2 c i 

Specializing to energy changes of the test particle due to 


elastic scattering, equation (40) gives 
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b l E 2 


, 1 n i n 2 6 f 

t * 61 2 K JU,')p°E° J [ 


da 


dv 8 v / d(cos 4 1 )(1+ cos ^ ) 

2 <u i' )p 2 E 2 j i 1 r r K ' dcosl1 ' 


X 


|_ n "< E 2 Pr +m 2 E i% /b . sl/2) ■ n+P 2 P i' E i,2 


(41) 


where Q ± = exp (-b , s 1 ^ y" ) ± exp (-b s 1 ^ y + ). The kinetic energy 
c 1 c 1 c 

density e for the test particle is just (E^-m^)/n^. 

A simpler expression for the rate of change of LS energy can 
be derived by considering the interactions between two isotropic 
relativistic MB gases at different temperatures. The instantaneous 
LS energy loss rate of one of the gases, due to scattering from 
particles in the other gas, is given through equations (10) and (37) as 


1 n l n 2 U i U 2 

’mb 2 (l+6 u )K 2 (n')K 2 ftlp 


J dv r lybl )f dE* j 

1 *'m. J 1 
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du e 


■ b i E r b 2 E 2 r 


(42) 


Y (E. -E ) 
c i 


vl- -U 


+ P V (p. - p CO s 4* ) u 

C C 1 


The Kronecker delta function accounts for the self- interacting case, 
and azimuthal symmetry is assumed in the scattering process. 
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Equation (42) can be simplified by recalling the definitions of 

and E^ from equations (32) and employing the following substitu- 

1/2 1/2 

tions: y =ch9.; (m y +m )/s = ch0_;(m 1 Y +m.)/s = 

c 1 2r 1 2 1 r 2 

ch9^ . It is then necessary to calculate a pair of two-fold integrals 

Li 

along the lines suggested for the calculation of the integral in equa- 
tion (35). These are 


'./■ 

1 


I j = / d9 j ch 9 ^ sh 2 9 , / du e e ^ U . and 
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(43a) 


" 1 

= f d9jSh^9^ f duue e^ U , 
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'-1 


(43b) 


where H = chQ ^ (u ' chO^ - u ^ ) an d ^ = shO^ (|a ' sh0 2 - |a^ sh0 2 )• 

1 /2 

Identifying -H ± 6 with -z ch(9 T B), one finds that shB = 

- 1/2 

z (Mj shG^ “ ^2 s ^®2^* Using the integral representation for the 

Bessel functions and the recursion relation K , (x) - K , (x) = 

v + 1 v- 1 

2v K (x)/x (Abramowitz and Stegun 1965), one obtains I, = 
v 1 

2chBK^(z ^ 2 )/z ' and I^ = -2shBK^(z ^ 2 )/z^ 2 . Equation (42) 
become s 
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The hyperbolic functions are given explicitly by 

1/2 

sh B = m 1 m_ (3 Y (b - b )/(zs) and 
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(45b) 


Though equation (44) is still rather complicated, it reduces in 
a number of cases. For example, if the scattering is elastic, 

rp »■(» 

d CT /dE d cos lb -» 6[E - E. (v )](da /d cos 4* ), and we recover 

l r 

the formula derived by Stepney (1983), namely 
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(46) 


If the gases are at the same temperature, then shB - 0 and chB - 1. 
Therefore 


n l n 2 M l ti 2 


/ 2 _. T _ * 1/2. „ 
f (y - 1 )K (bs ) r o./ * 
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(47) 


The second integral is simply e (y )ct(y ), defined in equation (27), 

.u 

* 0 ' T * * T * T 

because E - E = E^ , where E^ is the average energy of a 
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secondary created in an inelastic collision, so equation (47) properly 
reduces to the LS energy density loss rate -- the luminosity -- due 
to particle collisions of one isotropic MB gas with a second at the 
same temperature. The kinetic energy density of a relativistic MB 
gas is given by 


4tt A 


a 

/ 


dp p 2 (E 


-bE _ 

m)e - nm 


3 K l (|i) 

— + — 1 

M K„(M) 


(48) 


which, with equation (44), permits the calculations of the energy loss 
timescales. 

Calculations of the various timescales will be performed in 
the following chapter. Note that the energy loss timescales com- 
puted through equations (44) and (48) refer to the average energy loss 
for the entire distribution function whereas, in fact, one portion of 
the distribution may fall out of thermal form more rapidly than 
another. In this sense, Gould's treatment cited earlier is more 
accurate since it considers that portion of the distribution that would 
be first to fall out of equilibrium or last to be thermalized. This 
treatment does however have the advantage of applying to all temper- 
atures, and particularly the important transrelativistic regime. 


E. Compton Scattering 

The process known as Compton scattering, in which a photon 
is scattered by an electron, is of considerable astrophysical 
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importance. Felten and Morrison (1966) calculate the diffuse X- and 
y -ray background emission from the "inverse" Compton scattering 
of ambient low energy photons by the high energy cosmic ray elec- 
trons. Inverse Compton scattering is also proposed to explain the 
hard X-ray spectrum of Cygnus X-l in a two -temperature accretion 
disk model (Shapiro, Lightman, and Eardley 1976). Numerical 
studies of the Comptonization of low energy photons by a cloud of 
relativistic electrons have been used to model emission from X-ray 
sources (Pozdnyakov, Sobol', and Syunyaev 1 978; Sunyaev and 
Titarchuk 1980). 

Jones (1965) derives an expression for the rate of change of 
energy of an electron Compton scattered by a monoenergetic photon 
gas, and applies the results to the energy loss rate from a Planckian 
distribution of photons. The exact spectrum of Compton scattered 
photons for an isotropic monoenergetic photon source is also calcu- 
lated by Jones (1968); the complexity of the exact result renders it 
effectively useless except for the purpose of verifying limiting forms. 
Approximate expressions for the Compton power and the Compton 
scattered spectra in the Thomson and the extreme Klein-Nishina 
limit can be found in the review by B lumenthal and Gould (1970). 

The purpose of this section is to briefly treat the Compton 
power in the regime appropriate to the problem of Compton scatter- 


ing between a relativistic electron gas and an internally produced 
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photon source. The scattering oftentimes takes place intermediate 
to the Thomson and Klein-Nishina limits. Results will be applied to 
energy and momentum exchange rates between a relativistic photon 
and electron gas for the purpose, for example, of calculating the 
Compton- driven radiation forces which have been suggested as the 
mechanism for driving the jets observed in the vicinity of some 
galactic nuclei. Later uses of the formalism will be for the calcula- 
tion of the average energy increase of a photon per Compton scatter- 
ing due to an isotropic MB electron gas, in the manner of Lightman 
and Band (1981). An approximate Compton scattered spectrum can 
be calculated for such systems. 

The formalism used in Sections B and C to derive equation 
(26) for calculating production spectra from interacting gases of 
particles can be modified to a photon-electron system. The boost 
transformation is now to the rest frame of the electron rather than to 
the CM of the two particles. Letting y be the Lorentz factor of the 
electron in the LS and 9 the angle between the photon and electron 
momenta in the LS, define 

x = yk (1-8 cos 8 ) , (49) 

where k is the photon energy measured in units of the electron rest 
mass. The quantity x is the dimensionless photon energy in the 
CM. If x « 1, we are in the Thomson limit, whereas if x » 1, 
the scattering is in the Klein-Nishina limit. 
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From equation (21), the spectrum of Compton scattered pho- 
tons with energy E (also measured in units of the mass of the elec- . 
tron) due to a scattering characterized by the reaction parameters x, 
■y , and k is 

T 1 

d0 J d(cos 4 1 ) F^(E » cos ^ » x), (50) 

(cos 40 . 

mm 

v 'i' 

where E and cos 4* are given in terms of E, y , and cos 4 1 
through equations (18). The invariant Compton cross section is 
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(51) 


(Jauch and Rohrlich 1976), where r is the classical radius of the 

e 

electron. 

Following procedures similar to those leading to equation (26), 
the Compton scattered spectrum becomes, after integrating over 
cos 0 in the delta function. 
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defining z = 1/x + 1 - 1/E^, C = uv, and D = (1-u^) (1-v^) , 
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with u = (k - yx)/8yx and v = (E - yE )/QyE . integration of 
equation (52) is complicated by the appearance of E in v, but has 
been accomplished, with an additional integration over x from 
yk(l - p) to yk(l +3), corresponding to an isotropic monoenergetic 
photon distribution, in the paper by Jones (1968) just cited. 

Now define 


(E n a ) 


CO 



dE E n 


dcr(E; x, y, k) 
dE 


(53) 


as the n moment of the scattered spectrum. Thus, (E°a) = CT ^(x) 
is the exact Compton cross section, (E*a) / (E°a ) defines the 
average energy of Compton scattered photons in a scattering process 
characterized by x, y and k, and so on. Substituting dcr /dE 
from equation (52) into equation (53), and using the relation 
E = yE (1 + (3v), we have 
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dE (yE f 
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2 2 2 1/2 

where R = [(l-u - z ) + 2vuz - v ] . The limits are given by 

v^ = uz T (1 - u^)^^(l - z^)^ ^ from a consideration of the scattering 
process. Defining 


h 'i- 


dv v 
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(55) 
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straightforward integration gives 1^ - tt, 1^ = uzTT, I = 

2 2 2 2 

[3u z + (1-u -z )]rr/2 , etc. For n = 0, one obtains the Compton 
cross section 
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nr 

(E°ct) = Og (x) = — ~ 
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(56) 


reducing to the appropriate limiting values (Heitler 1954). The 

2 

Thomson cross section o = 8-rrr /3 . 

i e 

The first moment, the average scattered photon energy, is 
given through 

<E> = <E 1 a)/a c ( x )= [1 - T (x)] + kr (x) (57) 


using the relation k - yx(l +Bu). The function T(x) in equation (57) 
is defined as 
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and has the limiting forms 

! 1 - x , x « 1 , and 

4/3 2x , x » 1 . 


(5 9) 
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Since the initial photon energy is just k, the change in the photon 
energy due to Compton scattering is 

AE = ( E ) - k = [1 - T(x)] • (yx - kx - k)/x . (60) 


The integration over the variable x appropriate to an isotropic dis- 
tribution of photons in order to determine the exact electron energy 
loss rate can be evaluated (Jones 1965). For the present purposes, 
it will be more useful to consider approximations based on the be- 
havior of functions related to T(x). The appropriate functions, 
namely a c (x)/<J T , u ^(x) [1 - T (x)] /a ^ , and a c (x) • r(x)/a T are 
shown in Figure 2. 

The energy exchange rate due to Compton scattering between 
a photon and an electron gas is given by equations (10) and (60) as 


i /% h ( y / d 


V-el = t / dn ph ( ^ l dn el ( Z } I d(cOS 0) • AE • 0 >- V X) 


2 3 


CD 

f d ^ f el (Y) f 

yi 


= Sir m A ph A el / dk f ph (k) I dy f el (y) / dxo^x) 

J 0 


yk(l+ 8) 

yk(l- B) 


• [1 - T (x)] (yx - k - kx) , (61) 


where the last expression applies to isotropic photon and electron 

gases. Equation (61) properly reduces when the scattering is in the 

Thomson and Klein- Nishina limits (Blumenthal and Gould 1970). 

From Figure 2, one sees that for x between about 0.05 ahd 300, 

* 

o ^{x) [1 - T(x) ] /ct^, is equal to about 0. 07, within a factor of 2. 



Figure 2. Associated Compton power functions. The dimensionless 
functions a^(x)/aj (label 1), CT^(x)*r(x)/CT,j. (label 2), and 

a c (x) [l - T(x)] / Orjy (label 3) are graphed along with some 

approximating formulas. Functions are defined in the text. 
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/ 

Therefore the energy exchange rate is given approximately by 


Jt 

ph-el 


0. 07 



(p)( 8^y - k) 


= 0.07o t [n ph <0 2 v> -n el <k>] , 


(62) 


the last expression referring to an isotropic distribution of photons 
and electrons. Here (k) is the average photon energy of the photon 
gas and (p y) is the weighted value of 8 y over an isotropic elec- 
tron distribution. 

Equation (62) is adequate for calculating the energy exchange 
rate between a relativistic electron gas and a gas of photons with 
energy comparable to the rest mass energy of the electron. For 
example, given an electron MB plasma at about 0.5 MeV, the rela- 
tive number of electrons with y ^ 10 is insignificant (« 0. 35%). 

All photons with energies between about 25 keV and 8 MeV are cov- 
ered in this approximation, and the Compton energy exchange rate is 
then approximately given by equation (62). 

For x « 1, an expansion of equations (56) and (60) gives 

2 

x •(J^(x) • AE = a rp(~ 1 +2 lx/5) • [x (k - y ) + kx ] , A calculation of the 
energy exchange rate gives 


f = a 

ph-el T 


/ dn phfc>/ dn el ( E.> j k (t P 


2 2 
Y 


( 63 ) 


, 2f 16 fl 2.42 2 31 21 3 2 2 | 
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reducing in the case of Thomson scattering of low energy photons to 
the term linear in k (O'Dell 1981) and extending and correcting a 
higher order expansion (Calvani and Nobili 1982). To calculate the 
force on a relativistic plasma along the direction of the photon beam, 
a problem of interest in the study of relativistic jets (Cheng and 
O'Dell 1981), note that the momentum change along the direction of 
the photon momentum k is Ak = E cos y- k, where cos y = 
cost) 1 cos 9 + sin ^ sin 9 cos 0, with cos 9 = (yk - x)/ f3yk and 
cos 4 1 = (yE - E )/8yE, Calculating the rate of change of momentum 
in a manner similar to the earlier calculation for the rate of energy 
exchange gives 



1 +2x 


• [z sin 9 sin 9' + y cos 9 O + z cos 9 7 )] 



(64) 

and cos 9 = (k - yx)/Byx. 

Higher order moments can be calculated from equation (54), 
although the algebra becomes increasingly tedious. An expansion in 
moments can be used to calculate approximate Compton scattered 
spectra. In Chapter V, we will use the approximation based on the 
average value of the photon energy after Compton scattering by a 
relativistic MB plasma to treat Comptonization. 



III. APPLICATION TO PROCESSES IN A 
RELATIVISTIC PLASMA 


A. Introduction 

The principal use of this work is to calculate model spectra 
for astrophysical systems composed of relativistic plasmas. The 
most obvious example of such a system is a stable relativistic 
electron/ion gas. Because of pair and photon production mechanisms, 
a positron density will build to a steady- state value determined by the 
competition between pair production and pair annihilation. In addition 
to pair production from particle -particle collisions, photons will con- 
tribute to positron production through photon-particle and photon- 
photon pair creating processes. 

In the absence of a magnetic field or an external photon source. 


the ambient photon density is a function of the optical depth T 


(n, + n )cr (k) L ~ n <7 L where L is a characteristic linear 
+ - pn - 1 

dimension of the system, so we should guess that the steady state 


positron density n + is a unique function of the temperature and the 
optical depth of the gas. The optical depth T is technically a func- 


tion of the sum of the electron and positron densities, but since posi- 


trons are produced conjointly with electrons, it is sufficient to identi- 
fy (n + + n_) with n = n for order -of- magnitude estimates of the 
various processes, as shown in Table 1. The characteristic linear 


dimension of the system is denoted by L, and the sum of all photon 
absorption or scattering cross sections cr^Ck), a function of the 
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Table 1. Principal Processes in a Relativistic Gas 


Pair Producing 
Processes 

Optically Thin 
(T « 1) 

Optically Thick 
(T » 1) 

(Order of Process) 

+ - 

ee -» eee e 

2 -1 
n + ~ a n t c 

2 

a 

+ - 

ye -» ee e 

2 -1 
n + ~ a n T t c 

a 

+ - 

yy -» e e 

2 2-1 
n ~ a nT t 

"T* 

1 

Photon Producing 
Processes 



+ - 

e e -» yy 

n ph ~ n + T 

1 

ee -» eey 

n . ~ a nT 
ph 

a 

e y -* eyy 

2 2 

n ~ a nT 

ph 

2 

a 

Energy Changing 
Reactions 



ey ey 

2 

n , ~ anT 
ph 

1 

yy yy 

4 3 

n , ~ anT 
ph 

2 

a 

The r malization 
Process 



ee -» ee 

-- 

T 

Note: The symbol 

"e" stands for particle, 

"y" for photon, and 


tp, = na T c . The order of the process in the optically thick 

2 2 2 

case is compared to a ~ r = a , where X - _ is the . 

6 v> C> 

electron Compton wavelength. 
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photon energy k, is replaced by a value characteristic of the chief 
scattering process, namely the Thomson cross section for Compton 
scattering. 

The dominance of the various processes depend on whether the 
gas is optically thick or thin, as indicated in Table 1. In the Table, 
all processes with four vertices or less in the Feynman diagram are 
included, excluding associated radiative processes which are in all 

cases of order a (= 1/137) smaller than the basic process (for 

+ - + - 
example, y y -» e e y will always be dominated by y y e e , so 

maybe neglected). However, the "double" Compton process 

ey -» ey y is retained, as it can be an important source of low energy 

photons. 

In the extreme optically thick case, the plasma achieves its 
fully thermalized form with zero chemical potential for the photons, 
electrons, and positrons. In this event, the number densities of the 
three components are comparable in value, so that the dominant 
processes are determined by the lowest order Feynman diagram. 
Hence the most important processes are, beside the thermalization 
mechanisms, pair production through photon-photon interaction, 
electron-positron annihilation, and Compton scattering. 

In the optically thin case, it is necessary to determine the 
relative importance of the various photon and positron production 
processes. As T -» 0 , pair production through particle -particle 
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collisions is the dominant pair producing mechanism. The photon 

density n , ~ n,n_R (t) . , where R [cm"Vsecl~ a} o c is the 
1 ph 1 2 re ph ’ re J T 


reaction coefficient, i is the order of the process, and (t) L/c 

is the average photon residence time. Therefore the annihilation reac- 
+ - 


tion e e -» yy and bremsstrahlung production of photons are of pri- 
mary importance. If the optical depth is moderate (T ^ © (0. 1)), 
Compton scattering can also be significant, as can synchrotron emis- 
sion, if a sufficiently strong magnetic field is present. 

Thermalization processe s are primary in both cases as they 
determine the character of the distribution. If they dominate all in- 
elastic processes, the particles obtain a thermal form. The conditions 
under which a MB distribution can be assumed is treated in the last 
section of this chapter, and the remainder of the chapter is devoted to 
calculations of production rates and spectra on the basis of the formal- 
ism established in Chapter II, assuming that the particles achieve MB 
form. Also, opacities and absorption coefficients will be treated in 
Section D so that the proper photon density can be determined. 


B. Electron- Positron Annihilation 

The cross section for electron-positron annihilation in the 

frame of reference in which one of the particles is at rest and the 

other has Lorentz factor y can be written as 

' r 
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tt r 


C (y ) 
a r 


(y + 1) 

r 


y + 4y + 1 
r r 

v 2 -l 


My + (y 2 - D 1/2 ] - 

r r 


Y +3 

r 


(Y 2 - D 1/2 
r 


( 1 ) 


(Jauch and Rohrlich 1976). The classical radius of the electron is de- 
noted by r and the limiting forms of o (y ) are given by 
0 3. r 


2 

o (y ) > TT r / 6 

a r e r 


CT (y ) > TT r^ 

a r e 



Y - 1 « l ; 

r 

Y » 1 . 
r 


(2a) 


(2b) 


As measured in the LS, the total energy released from a pair of 
particles with LS energies E^ and E^ before annihilation is, from 
equation (II-6), 

e (Y r , Y c , u) = E 1 + E 2 = Y c m e [2 ( Y r + D]^ , (3) 

where m^ is the electron mass. If the gas is at a temperature great - 

6 

er than about 10 K, annihilation proceeds primarily through the emis- 
sion of two photons of equal energy in the CM of the annihilating parti- 
cles (Crannell et al. 1976). Thus the energy liberated per photon as 
observed in the CM is 

e (y ) = (m /2) [2 ( y + l)] 1 ^ 2 (4) 

r e r 


2 *2 2 . , * * * 

since s - 4m y - 2m (y +1) and E„ = m y , where y is the 
e e r 0 e ’ 1 

electron (or positron) Lorentz factor in the CM. 
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Noting that the positron and electron masses are equal, though 
the particles are different, the reaction rate (11-16) becomes 


r (|i) = 
a ^ 


C 

n n + U c C 

K 2 (M) A 


K ] [M[2( y +1)] 1/2 } 
dy (y - 1 ) rrjz a (y ) 

r r [ 2 ( v r ^)] 1/2 " 1 


(5) 


defining |J = m^/k^T and letting n ^ represent the electron (posi- 
tron) density. Using equation (3) and the reasoning associated with 
equation (II- 10), or equation (11-28) and (4), the luminosity is 




n n + n c 

K 2<H> 


a 

/ 


? 1/2 
dY (y - 1) K 0 f |4 [2 ( y + D] ' } o ( y ) 
r r 2 r J a r 


( 6 ) 


The result for the production spectrum follows from equations (11-26) 
and (4), namely 


S (k;M ) 
a 


n n, |ic 

- T 


- M k 


2k 


7 


dy (y -l)e 
r r 


M y r /2k 


ct (y ) 
a r 


(7) 


where k = E/m e , the dimensionless photon energy in units of the rest 
mass of the electron. (Equation [7] contains an additional factor of 
two because each annihilation yields two photons. ) 

The analytical expression (5) for the reaction rate is given in 
the work of Weaver (1976) and has been presented for the problem of 
the annihilation of electrons and positrons, together with the luminos- 
ity, in two papers by Svensson (1982, 1983). In the low temperature 
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(p >1) limit (though not so low that Coulomb corrections need be 

applied) and in the high temperature ( jji < 1 ) limit, agreement is found 

between approximations to the analytical results and the results of 

Zdziarski (1980) and the numerical work of Ramaty and Meszaros 

(1981) (though note a factor 2 error in Figures 2 and 3 of their paper). 

An examination of the spectral emissivity (7) reveals that in 

the low temperature limit, the position of the peak of the emission as 

2 

a function of temperature goes as m^(l + 3/4p + 9/32 p ), retaining 

1 /2 

the quadratic term. A second order expansion gives m^fe/p) 

(1 + 3/4p) for the full width at half maximum (FWHM) in this limit. 

At high temperatures, equation (7) agrees with the expression derived 
by Svens son (1983) using detailed balance arguments. The spectral 
emissivity approaches 

S a (k; U )-^n_n +w 4 kc ^ 2k ^ - (l-e« 2) e‘ u/2k 

( 8 ) 

for p « 1, where E^(y), the exponential integral function, is given 
by 

00 

E i (y ) = J dt e t t 1 *■ -y E - Bn y; y « 1 (9) 

y 

and = 0. 577... is Euler's constant. In the high temperature 

-EL/ 

limit, the energy of the emission peak is also given by Svens son 

2 1/2 

(1983) and the FWHM is given by 2 um^/p (u -u - 1) , where 
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u = 2 Bn {2 / jjl ) - y . 

The behavior of the FWHM derived through second order ex- 
pansions of the spectral emissivity about its peak reproduces the tem- 
perature dependence deduced in the numerical studies or in graphical 
representations of the approximations to the analytical results. It is, 
however, systematically 10-15% low in both limits, indicating that 
higher order terms are important. Thermal annihilation spectra at 
selected temperatures are shown in Figure 3. Note the broadening 
and blueshifting of the peak with increasing temperatures. 

Equation (7) is used to describe the annihilation spectrum of a 
relativistic MB gas composed of electrons and positrons. In the ex- 
treme optically thick case, the proper distribution function is instead 
a Fermi-Dirac distribution, so equation (7) will not be appropriate for 
such a system. Ramaty and Meszaros, in the paper previously cited, 
also numerically calculate the production spectrum for a Fermi- 
Dirac distribution of electrons and positrons, and find that the reac- 
tion rates and peak positions are within 10% of the values for a MB 
distribution at temperatures less than 10 K 0.85 MeV), with the 
greatest divergence of the spectral shape in the high frequency tail. 

C. Brems strahlung 

Bremsstrahlung is the principal photon producing mechanism 
in relativistic gases in the absence of a magnetic field. The object of 
this section is to present numerical calculations of the thermal 
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Photon Energy k 


Figure 3. Thermal annihilation spectra. The thermal annihilation 
spectra from an electron-positron plasma are shown for 
temperatures u~ l =k B T/m e = 0.1, 0.3, 1, 3, and 10. 
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bremsstrahlung spectra for MB electron-ion (e-i), electron- electron 
(e-e), and electron-positron (e -e ) gases. The positron-positron 
bremsstrahlung emission is equal to that of e-e bremsstrahlung, and 
the positron-ion bremsstrahlung is the same as e-i bremsstrahlung. 

Except for the fact that the differential cross sections can be 
complicated, the calculation of thermal bremsstrahlung is straight- 
forward. It is customary to present the results in terms of the ther- 
mal averaged Gaunt factor g(k; (i), where (-1 is the inverse tempera- 
ture m^/kgT and k, as before, is the dimensionless photon energy. 
The bremsstrahlung production spectrum S(k), equation (11-26), is 
related to the Gaunt factor g(k; [i) through the relation 


S(k) = 


2n l n 2 ao °T 

k 



exp(- (Jk) g(k; i_i ) . 


( 10 ) 


For the case of e-i bremsstrahlung, the relevant differential 
cross section is the Bethe-Heitler (B-H) formula (Bethe and Heitler 
1934; Jauch and Rohrlich 1976), which is the lowest order Born approx- 
imation result. Thermal e-i bremsstrahlung is of considerable astro- 
physical importance, and expressions for it at nonrelativistic temper- 
atures (in terms of the Gaunt factor) can be found in the paper by 
Karzas and Latter (1961). Gould (1980, 1981a) calculates the high 
temperature, though subrelativistic, thermal e-i bremsstrahlung and 
Quigg (1967, 1968) calculates the same for extreme relativistic MB 
gases, both employing approximate forms of the B-H formula. 
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We calculate the thermal e-i bremsstrahlung using the exact 
B-H formula. It is interesting to note how the general formula, equa- 
tion (11-26), reduces to a form that can be deduced through elementary- 
considerations. The e-i thermal bremsstrahlung spectral emissivity 
is 


S . (k ; p, ) 
e-i 


n n.n|i . c 
e x i 


2K (p)K (n.) 
2 2' i 


( 

/ 


dy (Y - 
r r 


x b[s(v r )] 


1) 

T/2 



e-i 


dk / 

k' 


1/2 

dCT B-H^ k/ ’ Y r^ -by c (k,k 7 ) [s (v r )] 

dl? * 6 


(ID 


1 /2 

making some obvious notational changes. In equation (11), s (y^) = 

, 2 , 2 ,1/2 , _ r , 4 . ^ n/o 1/2 

(m. + m + 2m. m y) m. + m y , ? , = [s- (m. + m ) l/2s 

i e xer i er’ e-x i e 

« m (y - 1), y = (k 7 /k +k/k / )/2, and dCT„ (k 7 ,y )/dk 7 is the B-H 

° .D - irl ^ 


e r 


formula (proportional to Z for an ion of charge Z). Now p = 

m./k_,T » 1 and p = m /k„T ~ ©(1), as the ion mass m. » m 
l B e B l e 

1/2 1/2 

and k-rjT. Thus K (p .) -» (Tr/2p.) exp(-p .) and bs (y ) — 

B 2 x x i r 

p . + p y ^ . Reversing the order of integration in equation (11) and 
noting, because of the large value of p , that the exponential factor 
contributes to the value of the integral only in the immediate neighbor- 
hood of k = k 7 , we have 


CO 

, n e°i W [ , ,2 ,, ‘ MY r dl, B-H (k > Y r ) 

■i (k;u)= ^oir / V Y r- lle — ak 

2 1+k 


( 12 ) 



Equation (12) can be immediately derived by assuming that the protons 
are at rest in the LS, as advertised. Defining 
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2 2 
ar z 
e 




da B-H (k> V r ) 
dk 


(13) 


equations (10) and (12) give 


«e-i (k ^ ) 


_3_ 

16 


3u 

2tt 


1/2 


Z 2 c 

k 2 (m) 




-|iY 


dy (y - 1 ) e 
r r 


B 




1 + k 


(14) 


Values of the Gaunt factor for thermal e-i bremsstrahlung are shown 
in Figure 4a at selected temperatures. 

The calculations of the production spectra for thermal e-e 
and e -e bremsstrahlung are much more difficult because of the 
extreme complexity of the cross sections. The e-e bremsstrahlung 
cross sections is given to lowest order by Haug (1975a) in terms of 
the cross section differential in the energy and the angle of the out- 
going photon. He also (Haug 1975b) calculates thermal e-e brems- 
strahlung at 100 keV, The cross section for electron-positron brems 
strahlung is evaluated by Swanson (1968), to lowest order. It re- 
quires two integrations over the angular distributions of the outgoing 
particles and an additional integration over the angular distribution of 
the produced photon to get it in a form suitable for equation (11-26). 

Defining, in accordance with equation (13), 
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Figure 


4. Gaunt factors for bremsstrahlung from relativistic thermal 

plasma. The Gaunt factors are given as a function of yk , 

where u * = 0 = k^T/m^ is the dimensionless temperature 

and k is the photon energy in units of the electron rest mass 

energy, for values of the parameter © = 0. 1, 0. 3, 1, 3, and 10. 

The Gaunt factors for thermal e-p bremsstrahlung, g , are 

e-p’ 

shown in Figure 4a, and the Gaunt factors for thermal e-e 
bremsstrahlung, g , are shown in Figure 4b by the solid 

G ~ 0 

lines. The values of g + /2, where g + is the Gaunt 


e-e e-e 

factor for thermal e + e” bremsstrahlung, are shown by the 

dashed lines in Figure 4b. At 0 s 1, g = 2g to the 

e -e 

limit of accuracy of the integrations in the approximation used 
(see Appendix), over the range depicted. 
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a r do (k. y ) 

e , /, x ee r 

— 1 (k, y ) - — 

k ee r dk 


(15) 


equations (11-26), (10), and (15) give 


/ \1 /2 


3k 


Uk 


32(1+6 12 } K 2 (M) 

2 1 
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!. 


dy 


(Y - 1) 
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r [2(v r + D] 1/2 
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f. 


ee „ / 


f (k' , y ) exp 
, /2 ee ’ r 


- M 


Y + 1 
r 


1/2 


Jl + Jl 
T k 7 


(16) 


where = (Y r ~ 1 ) / [2 (y + 1] and 6^ equals 1 for e-e brems- 

strahlung and 0 for e -e bremsstrahlung. 

Thermal e-e bremsstrahlung has recently been calculated 
(Stepney and Guilbert 1983), and agreement is found with this calcula- 
tion. The Gaunt factors are given, for various temperatures, by the 
solid line in Figure 4b. Thermal e -e bremsstrahlung has not been 
previously calculated. The calculation is performed in an approxi- 
mate fashion, as described in the Appendix. The Gaunt factors for 
this process are shown in Figure 4b by the dashed lines. 

In terms of the function f(k, Y ) defined for the three types of 
bremsstrahlung, the luminosity is given by equations (11-27) and 
(11-28) as 


i 


2 2 

n,mn, U_ ar Z c 



< 1+5 12 )K 2<WV 



(Y^-1)K 2 {b[s<Y r >] 1/2 
[b [s< Yr )] 1/2 } 



? 

dk 7 f(k' y ). 
» r 


( 17 ) 
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Alternately, it can be determined from an integration over the dimen- 
sionless energy k, weighted by k, in the production spectrum (10). 
Calculations of the luminosity are useful for the evaluation of energy 
loss timescales in Section F. 

D. Absorption and Pair Production 

In this section, the connection is drawn between the absorption 
coefficient and the reaction rate for a particle or photon traversing a 
relativistic particle /photon plasma. Application to sources of opacity 
such as Compton scattering and pair production are made. At suffi- 
ciently low photon energies, bremsstrahlung absorption can also be 
important, and is also considered. 

The treatment will be brief, as much of the information can be 
found in other papers dealing with the subject, for example, Lightman 
and Band (1981), Lightman (1982), Gould (1982a), Svensson (1982), 
Zdziarski (1982), and Stepney and Guilbert (1983). 

1. Opacity Considerations 

The absorption coefficient for a particle or photon of momen- 
tum p passing through an isotropic particle or photon gas described 

** w 

by the differential momentum distribution dm (p 1 ) is just the proba- 

1 MV 1 

bility for scattering (or absorption, as in a pair production process) 
per unit path length. Designate the absorption coefficient by K (p ). 


Then 
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(p 2 } = // dn i { Pi )da • V ( 1 - ( Sr£ 2 ) 


(18) 


This is related to the reaction rate r, equation (II- 1), through 


7- 


r = c( 1 + 6 12 ) /dn 2 (p 2 ) K^) . 


(19) 


Now 


introduce the invariant quantity q = ^ ^ = ^1^2 ” 


PjP 2 cos 9 , in analogy to the earlier use of y^ , defined in equation 

2 

(II-5). The utility of q over y^ is that problems associated with 
massless particles (i. e. , photons) are thereby avoided. For a par- 
ticle or photon of momentum p traversing an isotropic gas of pho- 
tons, electrons, positrons, and ions, each with distribution function 
described through equation (II-4), equation (18) becomes 


2tt A 


K (p ? ) = — f 
2 P 2 E 2 


r° r E l E 2 +p l p 2 

1 / J ,, 2, . 2.. 4 2 2.1/2 

-I dE i f i ( E i ) / d(q )a(q )(q -nijir^) 

i J TT» 


El E 2 -PlP 2 


( 20 ) 


Equation (20) can also be written as an integral over , 

+ 

2TT ^1 f 2 24 2 2 1/2 f ^ 

Mp 2 )= 7e / d(qV(q )(q -m^mp / dE i ^ 0^ )) , (21) 

2 2> t-, "■ 


m l m 2 


where 


_± _ 2 , 4 2 2,1/2 ., 2 

E 1 " [E 2 q * P 2 (q ' m l m 2 ) l /m 2 


( 22 ) 


As an illustration, consider the "absorption" coefficient in 
particle-particle scattering, when the particle of momentum p 
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traverses a relativistic MB distribution of particles described by- 
equations (II- 13) and (11-14). It is elementary to show that in this 
case, K (p ) = cr (p_), where r (p ) is the individual test particle 

La X, La t b 

reaction rate given in equation (11-33). 

Evaluating equation (21) in the limit m 2 ~* 0, that is, for a 
photon traver sing a MB gas, one obtains, defining k = p ? /m , 

La e 


K 

Y - par 


(k) 


n 


par 

2k 2 K 2 ( W ) 


dk 7 k 7 CT (k 7 ) exp 

y-par r 

thr 



(23) 


The quantity k 7 represents the photon energy in the rest frame of the 

particle, and k 7 is the threshold energy of the process. For the 
thr 

Compton scattering opacity, CT(k 7 ) -» a^(k 7 ) given by equation (11-56), 

and k 7 = 0. For low photon energies (k u « 1), the opacity is 
thr 


K c (k) = (n + + n_)a T 


k K_(u) 5 k K,(u) 


(24a) 


In the Klein- Nishina limit, that is, for k/u » 1, we find 


K c (k) 


(n, +n ) tt r 
+ e 

kK 2 (p) 



00 \ 

c ) 

K 1 (u)[«n(2k) +y] 

+ / dxxshx exp(-|jchx)> 


0 ' 


(24b) 


reducing to the result of Gould (1982a) in the high temperature 

(|J « 1) limit. The cross section for pair production through photon- 

particle reactions, which also contributes to the opacity, is given 


below. 
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For the process yy -» e e , the absorption coefficient for 
photon- photon scattering is 

2 ® 

TTr 

K (k) = (4TTA ) 

y-y k 2 y 

using equation (20). Equation (25) is written in terms of the function 
cp (k' • k) defined by Gould and Schreder (19&7). The number density 
of an isotropic photon source is, from equation (II-4) 


dk' f (k # ) cp (k' • k) , 



2„ Particle -Particle Pair Production 

The cross section for electron-positron pair production in 
particle collisions has been computed by the "equivalent- photon" 
method at extremely high energies (Bhabha 1935a, b; Williams 1935) 
and through the methods of quantum electrodynamics (Baier and 
Fadin 1971). The review of Budnev et al. (1975) summarizes the 
results of these and other calculations. For the process ee -» eee^e 
and ei -> eie + e , the pair production cross section is given by 

V<V = JF Z 1 Z 2 a r e r (2Y r* ‘ (is) (2Y r> + Bfe < 2 V r ) + cj , 

(27) 

for y ^ 100. The charges of the colliding particles are and . 

thr 

For the ee reaction, the threshold Lorentz factor y =7, B = -11, 

r 

and C~ 100. For the ep reaction, y = 3, B = 2. 6, and 

r 
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C ~ 40. Unfortunately there are no calculations or experimental 
data for pair production near threshold for these two reactions. 

The cross section for pair production in the collision of two 
heavy charged particles is reduced by a factor of the order of the 
ratio of the electron and proton masses, so has insignificant effects 
in relativistic electron plasmas. 


3. Photon- Particle Pair Production 

The cross section for the process ye -» ee + e , with thresh- 
old k^ r = 4, has been calculated by Haug (197 5a). He provides the 
following fit, accurate to better than 0. 3%, for the cross section 
(Haug 1981): 


o (k) 2 3 4 

— - = [5.6 +20.4(k-4) - 10. 9(k - 4) -3.6(k-4) +7.4(k-4) ] 

ar 


-3 2 

X 10 (k-4) , 


4 S k ^ 4.6; 


2 3 

0. 582814 - 0. 29842k + 0. 04354k - 0.0012977k , 


3. 1247 - 1. 3394k + 0. 14612k 
1 + 0. 4648k + 0. 016683k 2 


4.6 < k ^ 6 . 0; 
6 <: k <; 18 ; 


(28) 


28 „ 218 , 
= + 


1 T 4 3 2 

— - -&n 2k + 3. 86 3 On 2k - 1 1 Bn 2k + 27 . 9 , 

k L 3 J 


k £ 14. 


A fit to the cross section for pair production in photon-proton 
scattering is given by Stepney and Guilbert (1983), accurate to 0. 1%. 
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2 1/2 

Defining e = (k -4) /k, the fit is: 


CT (k) 6 2468 

— ^ — = a-^r- (1 +0. 87 5e +0.7 55 e +0. 661 e +0. 589 e ) . 
32 


2$k^ 2.4; 


= (7.620 - 8. 0218k + 2. 520k 2 - 0. 2047k 3 ) X 10" 4 


2. 4 ^k £ 4. 0 ; 


(29) 


(&n 2k - 109/42) + (473.65 +241. 26 Bn 2k - 81. 151 Bn 2 2k 
6tr 


+ 5. 3814 ^/i 3 2k) X 10~ 5 /k , k s 4. 0 


The threshold for the process is 2. 

The rate of pair production due to photon-particle collisions 
in a relativistic MB gas can be found from equations (19) and (23). 
The rate is 


r y-par 


4tt A 

Y 



f (k) K 
Y 


Y-Par 


(k) 


(Gould 1971; Weaver 1976). 


(30) 


4. Photon- Photon Pair Production 

The absorption coefficient for this process is given by equa- 
tion (25). From equation (19)# the pair production rate is 


16tt A A 

Y 1 Y 2 


r = (crrr ) - 

y — y e 1+6 


dk f (k) 


12 


03 

/ dk' f (k' 

/ Y 1 

•'l/k 1 


) cp (k • k' ) . 


( 31 ) 
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The Kronecker delta function corrects for double counting if self- 
interacting photon distributions are being considered. 


5. Bremsstrahlung Absorption 

Both bremsstrahlung (free-free) absorption and synchrotron 
self-absorption can be effective sources of photon opacity at suffi- 
ciently low photon energies and, for the latter case, high magnetic 
fields. From Kirchhoff's law, the bremsstrahlung absorption coeffi- 
cient K^k; I-0 f° r a photon of energy k traversing an isotropic 
relativistic MB gas is 


2 3 

n X C U k 

= T- <e MK - 1) S(k) , 

ck 


(32) 


where is the electron Compton wavelength ft/m^c (r^ = aX^.). 
Equation (32) can be written in terms of the Gaunt factor through 
equation (10). Because of the slowly- varying logarithmic depend- 
ence of the Gaunt factor at low frequencies, the principle behavior 
-2 

of goes as k as k -» 0. 

Inserting numerical values for K , it is readily seen that 

b 

for an optically thin relativistic MB gas, bremsstrahlung absorption 

is only effective for extremely low photon energies (k ^ &[l0 
10 3 

for n p~ I® cm )(Lightman and Band 1981). The synchrotron 

self-absorption coefficient K (k; \l) for an extreme relativistic 

s 

magnetoactive MB plasma has been treated by Jones and Hardee 
(197 9). though note Gould's (1982a) correction. 
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E. Particle Production in Relativistic Nucleonic Plasmas 

The inner region of an accretion disk surrounding a black 

hole has been proposed as a site of relativistic nucleonic plasmas, 

with temperatures in the range 0.01 ^ 0 = k T/m ^ 1 (Lightman, 

p B p 

Shapiro, and Rees 197 5). At these temperatures, pion- producing 
reactions from nucleon- nucleon collisions would be a source of 
energy loss and luminosity. Elements of such a system could only 
have a transient existence; otherwise pair production processes 
would conspire to transform the kinetic energy of the relativistic 
nuclei into an optically thick plasma of electron-positron pairs and 
photons (Bisnovatyi-Kogan, Zel'dovich, and Sunyaev 1971). 

The object of this section is to compile the experimental data 
required to calculate the magnitude of secondary particle production 
in accretion phenomena. We consider only plasmas composed of 
protons. The effects of heavier nuclei can be accounted for in an 
approximate way by employing a multiplicative correction factor 
(Orth and Buffington 1976). 

One might think that the cross section and production data 
would be readily available from cosmic ray studies. In fact most 
calculations of secondary production by cosmic rays in the past 10 
years have depended upon high energy data, where scaling laws are 
supposed to apply. It is the cross section values near threshold and 
at intermediate energies (laboratory proton kinetic energies E » 

K 
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1-10 GeV) that are most important for secondary particle production 
in accretion disks. Actually, these are also the most important 
energies for the production of the bulk of the secondary particles in 
cosmic ray interactions in the galaxy, but the scarcity of data and 
the promise of simple functional forms for the production cross sec- 
tions at ISR energies have led workers to give more attention to the 

❖ 

highest energy data. 

In addition, contrary to published assertions (Kolykhalov and 
Syunyaev 1979; Stephens and Badhwar (1981), some information on the 
momentum and energy distributions of secondary pions produced in 
proton-proton collisions at intermediate energies does exist. In 
Table 2 is tabulated all data between threshold and ISR energies 
which are concluded to have value as regards the production cross 
sections or the CM momentum or energy distributions of the pro- 
duced pions in proton-proton collisions. Numerical fits are evalu- 
ated for the pion- producing processes and are given below. Second- 
ary K mesons, and A and E baryons are neglected in this study, as 
they contribute at most 10% (by particle number) due to the small 
relative cross sections and large threshold energies. 

Distinguish four channels of pion production: pp -» TT + d , 
pp -» tt X, pp -» tt°X, and pp -> n X. The X refers to inclusive 


For example, Stephens and Badhwar (1981) show that 50% of all 
Tr°'s produced in the galaxy come from cosmic rays with E, ^ 
2. 5 GeV. " ' * 
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production (Section II. C) except in the case pp -» tt X , where it is 
understood that the contribution from the channel pp -a TT + d is ex- 
cluded, d standing for a deuteron. For energies near threshold, 
the cross section for pion production is described by the Mandelstam 
model (for a general reference, see Lock and Measday 197 0). At 
the highest energies, the inclusive cross section, denoted by cr 

TT X. 

for the process pp -» ttX, can be shown to have the form = 

A + C to s as s 00 , from elementary considerations of the scaling 

hypothesis. In this expression, s, defined earlier, is the sum of 

the four -momenta of the two colliding particles: s = 2m (E, + 2m ), 

p k p 

where E is the kinetic energy of one proton when the other is at 
K 

rest. Based on considerations of the fragmentation model of the 

production process, it is proposed in reference 25 to add the term 
- 1/2 

Bs . Since s ~ 2m p. at these energies, where p is the 

pi & 

momentum of one of the colliding particles in the frame of reference 

in which the other is at rest, the high energy functional form is 

- 1/2 

taken to be a^^fmb) = A +■ Bp^ + C 2n p^ . The coefficients A, 

B, and C are determined by minimizing the vertical variance of 
the fit with respect to the experimental cross sections. All mo- 
menta and energies have units of GeV, and the fits intermediate to 
the threshold and high energy regimes have no theoretical basis. 

Threshold laboratory kinetic energies for the four processes 

TT + d tt^X tt°X rr’x 

are E. = 0.288, E. = 0.292, E, = 0.280, and E, , = 

k, thr k, thr k, thr k, thr 
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0.600. The respective values of m^. (Section II. C) are m^ , 

m + m , 2m , and 2m + m ■ . The quantity T) is defined as 
P n p p tt + 


T] = 


' max 


[0 


2 — 2 

m - m ) 
rr x 


. 2 — 2 , 1/2 
4m m J 
TT x 


m_ 


rr 


2 m n s 


1/2 


(34) 


where m is the mass of the produced pion for the channel under 

consideration (m , = 0. 1396, m „ = 0. 135), and p is the maxi- 
ng tt° max 

mum pion CM momentum. Then the following fits are valid from 
threshold to the highest accelerator energies, within the accuracy 
of the experimental data: 


0. 18 ti + 0. 95 ri 3 - 0. 0 16 ti 9 , ^ 0. 65 


a (mb) = 
TT d 


-3 9 

0. 56 E, ^ , 
k 


0.65 £E k < 1.43 (35) 


0. 34 E. 


■2. 5 


1.43 ^ E, ; 

k 


0. 95ri 4 + 0. 099 t 1 6 + 0. 204 9 8 , p^ ^P 1 <0. 954 


c (mb) = 
TT X 


0. 67 ri 4,7 + 0. 3 , 


22.0(p -1. 27) 


0. 15 


0.954 < p < 1. 29 


1.29 < p < 2.81 


(36) 


•40 


- 1/2 


.9 + 57. 9p L + 27 . 0 fei p , 2.81 ; 



68 


a (mb) 
TT°X 


O (mb) 
tt _ X 


0/0 

0.032 r] + 0.040 1] + 0.047 r\ , 

32.6( Pl - 0.8) 3, 21 , 

5. 4(p^ - 0. 8)°* 81 , 

- 1/2 

-59. 5 +48. 5 +32.O0np 1 , 

2. 33(p - 1.65) 1 * 2 , 

0. 32 p 2,1 , 

- 1/2 

-69- 3 + 74. 2 p 11 +28. 2 Bn p , 


IT X 

Pl.thr ^ P 1 

0.954 ^p 


< 0. 954 

< 1. 29 


(37) 

1.29 £ p < 5.32 
5. 52 £ p ; 


1 .65 ^ p < 2. 81 
2.81 < p < 5. 52 (38) 

5. 52 £ p . 


The experimental data for the latter three channels of pion produc- 
tion, together -with the fits, are shown in Figure 5. 

From the fits (35) - (38), we can calculate the reaction coeffi- 

2 

cient R = r/n , where n is the proton density and r is the rela- 
P P 

tivistic reaction rate (II- 16). The results are shown in Figure 6 as 

a function of the dimensionless temperature © . For 0 « 1, the 

P P 

reaction coefficient is a factor of two smaller than found in two pre- 
vious calculations (Weaver 1976; Kolykhalov and Syunyaev 1979). 

The correctness of the present calculation was verified by using the 
cross section for tt° production given in the latter paper, confirm- 
ing their factor -of-two error. 

To evaluate the luminosity and production spectrum from a 
relativistic MB plasma of protons, it is only necessary to have the 



Table 2. Compilation of Pion Production Cross Sections in Proton- Proton Collisions 


0. 800-1. 30 

0.897-1.01 

1.06-1.29 1. 17 

1.16 

1. 29 

1.29 

1. 28 

1. 38 

* 

* 

0.91 ±0.15 

1. 2 * 0. 3 

3. 4 * 0. 4 


3.0 ± 0. 3 

3. 46 ± 0. 25 


♦ 

* 5. 21 * 0, 44 

5.0* 1.0 

10. 2*1.2 

13. 1 * 1. 2 

10. 8 t 0.5 



3.7 * 0. 3 4.3 * 0.1 7.74 * 0. 4 

18. 3 * 0.7 16.4 * 1.5 22. 2 *0.8 

0.01 ±0.01 1.27 5 * 0. 2 


2 


3 ,f 


2.8 
9. 35 ±0.4 
23. 4 * 0. 6 
2. 99 * 0. 2 


4.64 * 0. 2 


19.2*1 35.2*2.4 

29. 9 * 1. 5 42.7 ±0.7 

11.5*1 21. 1 * 0. 4 


8 ' 2 


24.0 69 

31.6 * 2. 6 41.6 * 3 53. 5 * 3. 1 

47.5 * 3 56.8 * 0.9 78. 5 * 1.3 

29.7 * 3 33. 8 * 0. 6 57. 5 * Q. 6 


80. 5 * 3 85 * 8 


91-8 ± 9. 2 108 * 11 125 * 13 

66.9 * 1. 3 86 * 2 99. 5 * 3 



22** 

22** 

22** 

22“ 

23 

Zi 


2 


107 * 7 


139 * 11 


158 * 24 169 * 24 


159.3*16.8 


Reference gives several cross section values. 

'Contains CM energy or momentum distribution data of produced plons. 
t Values given Include renormalization of cross section from Ref. 15. 


Values taken from compilation of cross sections. 
'Assessed errors for cross sections: *10%. 
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Figure 5. Inclusive pion production cross sections in proton- proton 
collisions. The experimental data of Table 2 for the inclu- 
sive cross sections, in millibarns, for the production of pions 
in proton-proton collisions are shown with numerical fits to 
the data given in the text. The data is plotted as a function of 
the momentum of one of the photons, p^ , when the other is at 
rest in the LS, Figure 5a: inclusive cross section for the 

production of TT°-mesons. Figure 5b: inclusive cross section 

for the production of tt + - mesons, excluding the channel 
pp h> TT + d. Figure 5c: inclusive cross section for the produc- 


tion of tt -mesons. 
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Figure 6. Reaction coefficient vs. © in a relativistic MB plasma 

P 

of protons. The reaction coefficient R for the production of 

tt - mesons in a relativistic MB plasma of protons is shown 

as a function of temperature ® = k_T/m . The reaction 

p B p 

coefficients for the channels pp -» rr^d (solid light line), 
pp -> tt°X (dotted line), pp -» rr + X (dashed line), and 
pp -» rr X (dot-dashed line) are shown separately, along 
with the total reaction coefficient (solid heavy line) for all 


channels of pion production. 



R [cm /sec 
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angle -averaged energy spectra of the produced pions in the CM. For 
the process pp -» rr + d, this information is automatically determined 
from two-body kinematics. Essentially complete information for the 
CM spectra of rr + and tt° are given in Refs. 4, 10, and 13 of 
Table 2. Reference 13 also gives the CM energy distribution of tt 
T o determine the luminosity, only the cross section weighted by the 
energy in particles produced, equation (II- 27 ), is needed. These 
numbers were calculated from the previously mentioned references. 
At higher energies, the values were derived from expressions for the 
Lorentz invariant production cross sections used in cosmic ray 
studies (Bahdwar and Stephens 1977; Bahdwar, Stephens and Golden 
1977). The average energy (E) for the three types of produced 
pions can be fitted to the expression 

<E> " * K(E k- E k ,thr> a ■ (3 ’> 

where K = 0 . 27 6 , 0 . 28 5, and 0. 27 5, and a = 0. 22, 0. 185, and 

0. 187, for tt , tt° , and n , respectively. The luminosity coeffi- 
2 

cients L = f /n for the various channels and the total luminosity 
P 

coefficient, calculated from equation (11-28), are shown in Figure 7 
vs. © . Also shown is the energy loss timescale, used to determine 
if a MB distribution of particle energies can be established in a rela- 
tivistic plasma of protons. Its significance will be discussed in the 


next section. 
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Figure 7. Luminosity coefficient and energy loss timescale vs. 0 

P 

in a relativistic MB plasma of protons. The luminosity 
coefficient L for the production of tt -mesons in a relativ- 
istic MB plasma of protons is shown as a function of temper- 
ature 0 = k_T/m . Key for the graph is the same as in 

P B p 

Figure 6, except that the curve labelled "T " refers to the 

hj i-j 

energy loss timescale due to pion production. The values on 
the right side refer to the energy loss timescale in units of 
seconds, and the values on the left side refer to the luminos- 

3 

ity coefficient in units of GeV-cm /sec. 
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F. Thermalization, Reaction, and Energy Loss Timescales 

We now treat the problem of the establishment of a MB dis- 
tribution in particle energies using the results of Section II. D. Con- 
sider first the case of a highly relativistic electron of Lorentz factor 
traversing a MB electron gas at temperature (in units of the 
electron mass energy) ® = [i ^ . Assuming » ®, the test elec- 
tron will be in the high energy portion of the distribution, which is 
the last to thermalize because of the reduction of the value of the 
elastic scattering cross section at high energies. The test particle 
will be thermalized if the characteristic time for elastic scattering 
is shorter than the timescales for energy loss or annihilation. 

The timescale for annihilation can be found from equations 
(11-33) and (2b). If the gas is subrelativistic, the annihilation time 
approaches 


a ’ ^ ^ 1 n (cTTr ^) [S/n (2y ) + 1 /2] 

- e c. 


(40) 


since y 0 »l, and y y (1 +j3 _ 8 ) -» 2y y and y y (1 - |3 (3 ) -» 
c. c r c. r c r c. r z r 

1 / Y r + in the limit "© « 1 inequation (11-33). For an 

extreme relativistic gas, the annihilation time is, in the approxima- 
tion y^ » ® » 1 , 

2y 2 © 

(cnr. e 2 ) [Bn (2y 2 ) + 1/2] 


T > 

a , t ® » . 1 


n 


( 41 ) 
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Compare this last result with the annihilation time computed 
for a MB gas composed of electrons and positrons at extreme rela- 
tivistic temperatures. Defining T =n,/r , with r given by 

r a, MB + a a 

equation (5), we obtain 


T > 

a, MB 0 » 1 


20 


n (cTT r ) Bn .9 
e 


(42) 


(The approximation requires the information that J dx 

00 . 1 
f 3/2 

exp(-x) = 1.13 and / dx fe x • x • exp(-x) = 1. 03. ) 


3/2 

x 

Comparing 


1 

equation (41) with (42), we see that (for this case at least) it is per- 
missible to identify y^ with 0 to get reliable estimates for the 
timescales from the test particle approximations. (Note the order of 
magnitude error for the result (41) cited by Lightman (1983). ) 

Thermalization in a pure electron gas occurs through Miller 
scattering. The cross section for this process is given by Jauch and 
Robrlich (1976). In the test particle approximation, the thermaliza- 
tion timescale can be evaluated through equation (11-41). Assuming 
» ® » 1* one can approximate « Q , y^ » 1, and retain 
only the leading terms in the Miller scattering cross section to get 


t, M^ 0 » 1 


n 


, 2 
(cfT r ) 
e 


^ 2 0 

f 1 - 

\8 


0n 2 - v)j 


m 


(43) 


In this expression, the minimum scattering angle is the ratio 

between the maximum and minimum momentum transfers, di = 

m 
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-17 -31/2 2 2 

hcu H ~xn ~ 4 X 10 [nn(cm )] /y, , where ou = 4rrr (ic n /3, 

p 2 e 2 p e e 

the square of the plasma frequency (Gould 1981). 

A similar calculation can be performed for a positron slowing 
in a pure electron gas through Bhabha (e^-e elastic) scattering. The 
characteristic time in this case is 


y 2 0 


t, Bha 0 » 1 2 

n (cTT r ) 2 

e 


77 - 071 ^ 
12 m 


(44) 


Equations (43) and (44) may be contrasted with the results of Gould 

*u 

••I s 

(1982a, b) by the inclusion of the Coulomb logarithm, -bn\b . Identi- 

m 

fying with |a ^ in these equations according to the argument fol- 

'I s 

lowing equation (42), and assigning for -bn 4^ the canonical value 

30 for astrophysical problems, it can be seen that both annihilation 

and thermalization timescales possess the same dominant behavior, 

2 

proportional to 0 . Considering the numerical coefficients, one 

finds that the thermalization time is always considerably shorter, by 
an order of magnitude, than the annihilation time, for © » 1. Thus, 
a positron should be able to thermalize before annihilating in an ex- 
treme relativistic MB gas of electrons. 

Whether the electron gas obtains the MB form depends on the 
timescales for energy loss of the electron gas. For the extreme 
relativistic case, Alexanian (1968) has calculated the electron- 
electron bremsstrahlung luminosity l . It is 
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2 2 T 51 

S. r — — * 24 n m ar c© 2n(2@) - y + — . (45) 

e-e © » 1 - e e L T E 4J 

For © >> 1, the kinetic energy density e approaches 3m n© . 

Comparing the energy loss timescale T = e /f with the 

6 ~ G K G ”6 

thermalization time T- . , equation (43), we find that for tempera- 

t, M0 

tures greater than about 10 MeV, energy loss through bremsstrahlung 
emission dominates the thermalization process, so the establishment 
of a MB distribution cannot be assumed. 

The thermalization of a gas containing electrons and ions 
depends on the energy exchange rate between the electron and ion 
component of the gas. Gould (1981) considers the energy exchange 
between electrons and ions through Mott scattering. For a neutral 
electron- proton extreme relativistic gas, he obtains 


T 

e-p 


3@_ 

© » 1 4tt 


m 


(n cr ) 
e 


— -—E- (fa 2 + 1 - Y - Bn vjj_) 1 
in m 


(46) 


where the temperatures of the electron and ion components of the gas 
are assumed to be of the same order of magnitude. Quigg (1968) de- 
rives the ©» 1 electron-proton brems strahlung luminosity 


SL 

e-p 


> 

©» 1 


12n n m ar c© U/i 
-pee 


© ^ 



(47) 


Note that it is almost exactly one-half the electron-electron brems- 
strahlung luminosity, equation (45), in this limit. The total energy 
radiated in an extreme relativistic e-p gas will be the sum of 
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equations (45) and (47). Proceeding as for the pure electron gas 
case, it can be determined that the characteristic time for energy 
loss by brems strahlung is always less than the thermalization time, 
equation (46). Thus an extreme relativistic MB electron- proton gas 
cannot be assumed to exist in nature. The timescales become equal 
at 1, that is, near 0. 5 MeV, but the relativistic approximations 

are no longer valid in this regime. However, numerical evaluation 
of equation (11-46), generally valid, yields a similar temperature 
(Stepney 1983). 

Elastic scattering between protons, including Coulomb and 
nuclear effects, has been treated by Gould (1982c) and Stepney (1983). 
The timescale for energy loss must be compared with the thermaliza- 
tion time by elastic scattering to determine if a MB distribution in 
particle energies for a relativistic plasma of protons can be estab- 
lished. The energy loss timescale is plotted in Figure 7. It was 
evaluated according to the relation 


T 


-1 

EL 


n (R 

P TT 



+ L n /e ) 
TT p k 


(48) 


where R is the temperature -dependent reaction coefficient for 

tt d 

the process pp -> tt d (Fig. 6), and L is the total luminosity coeffi- 
cient in pions excepting the pp -» TT + d channel. The kinetic energy 
density e of the relativistic proton plasma has been given in equa- 

K 


tion (11-48). 



To the accuracy of the experimental data (10-20%) and to the 
extent that pion production is the dominant energy loss mechanism, 
the calculated energy loss time is exact. Comparing the inelastic 
timescale to the elastic scattering time calculated in the previously 
cited papers, we find equality of timescales at 60 MeV. Thus, for 
temperatures less than 60 MeV, the establishment of a relativistic 
MB plasma of protons can in principle take place. 

If there were a magnetic field present, synchrotron energy 
losses could prevent the formation of a MB distribution. The ratio 
of the electron thermalization time , equation (43), and the character- 
istic synchrotron energy loss timescale for a relativistic MB dis- 
tribution of electrons (Gould 1982a), is 


T 

M^ > 

T © » 1 
syn 


128© 4 

3 



2 




(49) 


where e^ = B /8tt is the magnetic field energy density. Therefore 
at equipartition field strengths, the synchrotron cooling is able to 
dominate thermalization mechanisms and prevent the formation of a 
MB distribution. 

The ratio of the characteristic time for annihilation, equa- 
tion (42), and the energy loss timescale for synchrotron cooling, as 


computed above, gives 
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for 0 » 1. Hence a magnetic field two to three orders of magnitude 
less than the equipartition value can cool electrons from extreme 
relativistic temperatures before annihilation. This ratio is important 
for interpreting annihilation line features in gamma ray burst spectra. 



IV. PHYSICAL APPLICATIONS 


A . Boltzmann and Pair Balance Equations 

In this section, we examine the properties of static, non- 
magnetic, relativistic thermal plasmas. The appearance of such sys- 
tems in nature could conceivably come about in several ways: through 
dissipation of energy in shocks moving at relativistic speeds; mech- 
anical heating of a neutron star crust due to, for example, a phase 
transition in the interior (Ramaty, Lingenfelter , and Bussard 1981); 
the heating of a diffusive coronal plasma by an intense source of 
energetic photons (Guilbert, Fabian, and Rees 1983); accretion flow 
(either disk or spherical) near a collapsed object; etc. In what fol- 
lows, a number of approximations will be made in order to model an 
ideal type of system, namely a steady spherical plasma cloud of 
radius L with uniform density and temperature. This is the sort of 
object that Lightman and Band (1981) and Lightman (1982) have con- 
sidered at extreme relativistic temperatures 0 » 1 , and Svensson 
(1982) for all temperatures, but at small optical depth. In a recent 
preprint, Svensson (1984) examines the properties of this system for 
all optical depths and throughout the entire relativistic range of tem- 
peratures. 

Even the strongest gravitational fields, as on the surface of a 
neutron star, are inadequate to confine such a system (Colgate and 
Petschek 1981). A sufficiently strong magnetic field can, but the 
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production of photons through the synchrotron process changes the 
properties of the system, a circumstance beyond the scope of this 
work. In the absence of a confinement mechanism, the relativistic 
plasma will expand adiabatically, and a uniform density and tempera- 
ture cannot be assumed. Thus, the adopted model must be considered 
artificial, serving only to suggest properties of a more realistic and 
more complicated system. (Discussion of plasma confinement and 
trapping are found in Lightman (1982) and Svensson (1984).) 

Nevertheless, the assumption of homogeneity throughout a 
spherical region permits a quantitative analysis. Beyond this, we 
assume that the particles are described by a thermal distribution, 
which, as we showed in Section III. F, is not strictly true for © » 1. 
Qualitative behavior of solutions should however be expected to be 
correct, as only the magnitude of particle production and annihilation 
rates are affected, not the character of these rates. Also, photons 
that Compton scatter with energies k~ ©(©) are assumed to obtain 
a Wien form, as we examine in more detail below. External photon 
sources are not considered, and double Compton scattering is treated 
in an approximate fashion, as this process can be important for the 
calculation of the steady state positron density in the low tempera- 
ture, large optical depth regime (Svensson 1984). 

In general, Svensson's (1982, 1984) notation is used.' The 
density of protons is denoted by n^ , of positrons by n + = zn^ , and 
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two optical depths are distinguished, the proton Thomson optical 
depth Tp = n^a^,L, and the Compton optical depth t^Ck; 8 ) = 

K (k; 9)L , where K (k; 9) is the frequency-averaged Compton 
opacity, Eqs. (Ill- 23) and (III- 24), and is shown in Figure 8a for 
selected temperatures. All mechanisms associated with scatterings, 
absorption, or loss from the system are treated on the same footing, 
in terms of an associated absorption coefficient, or opacity. Thus 
we can write a differential equation for the photon number spectrum, 
n (k, t), namely 


3n (k) 

— — = -c { K (k) + K (k) + tC (k) + K [k, n (k)] 

3t L, C y-par y-y y 


+ K (k)} n (k) + S (k) + Q r (k,k') . 

b y y O 


( 1 ) 


S (k) is the photon source term associated with bremsstrahlung 

y 

processes, electron- positron annihilation, and double Compton scat- 
tering, and Q^fk, k 7 ) represents a photon source term due to 
Compton scattering of a photon from energy k 7 to energy k. The 
absorption coefficients in the brackets in Equation (1) are associated 
with loss from the system, Compton scattering, photon- particle and 
photon- photon pair production, and bremsstrahlung absorption, 
respectively. The temporal and temperature variables, as well as 
dependences on various system parameters, are implicit in the 
terms of the equation. 


\ 
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The form of K (k) can be established by examining the be- 

havior in the optically thick and thin limits. As (and T^,) -» 0, 

the average photon path length is about L; therefore K^k) 1/L in 

2 2-1 

this limit. For large optical depth, L = N^ m £p » where = 

K (k; 0) and N is the average number of "steps" required for a 

c 

2 

photon to go a distance L. Thus N = t^. and the average total dis- 
tance travelled by a photon is given by = ? T^L ' w ^ ere ? 

is a geometrical factor equal to about 1/3 for diffusive escape from a 
uniform spherical distribution of sources (Sunyaev and Titarchuk 

1980). We take 1C T (k) = n ct„/t (1+t„/3), bridging the two limits. 

Li pip L 

The other absorption coefficients are given in Chapter III, but note 

2 

that because the bremsstrahlung absorption coefficient n^ , 

whereas the rest of the coefficients are proportional to n , the num- 

P 

ber of variables required to completely specify the system are 3: 

0 , L (or Tp), and n^ . The steady state positron density is then 

z = z(©,t , n ) -* z (© , t ) for T ^ 1, because only when Compton- 
P P P P 

ization of extremely low energy photons is important should the effects 
of free -free absorption be felt on the rate of pair production. 

Previous studies of relativistic gases of photon and particles 
have proceeded by introducing characteristic photon energies which 
define regimes where the free -free absorption coefficient domi- 
nates K , and where it dominates K . A third photon energy is 
introduced which characterizes those photons suffering significant 
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changes in energy due to Compton scattering; i.e., where the so- 
called Compton y parameter exceeds unity (Shapiro, Lightman, and 
Eardley 197 6; Rybicki and Lightman 197 9). Semi -quantitative spectra 
are determined for the various photon energy regimes on the basis of 
the relative importance of scattering and absorption. 

Instead, we solve equation (1) directly by assigning to 
Q (k, k') an approximate form. Even so, it will prove easier to 
determine the spectra and properties of the system in various temper- 
ature and optical depth regimes by considering dominant mechanisms 
and solving equation (1) in an approximate fashion. 

Recall from Chapter II. E the expression derived for the aver- 
age energy change of a photon that is Compton scattered by an elec- 
tron (or positron) in a reaction characterized by the reaction parame- 
ters y» x » k. Averaging over a MB distribution of electrons, the 
exact average energy increase per scattering by a MB distribution of 
electrons and positrons is, from equations (11-56,57,58), 

00 

f dx e ^ a_(x)[r(x) - l] • [2|jk x +|jk - ux - 2xk] 

0 c 

A(k; 0) - 1 = , (2) 

2|jk f dx x a „(x) e ^ 

0 c 

/ x k \ 

where 0 = ~Z \"k + ~xj ' ^& ure 8b, function A(k; © ) is plotted 
for (i = 0 = 0. 1, 0. 3, 1, 3, 10. A photon with energy k before 
scattering emerges with energy A(k)*k after scattering. The exact 
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Figure 8. Compton scattering opacity and average Compton 

energy increase. In Figure 8a is shown the opacity of a 
photon with energy k due to Compton scattering by a MB 
distribution of electrons at various temperatures in the 
transrelativistic regime. The average energy increase 
of a photon that initially had energy k, due to Compton 
scattering by a MB distribution of electrons at different 
temperatures, is shown in Figure 8b. 
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average increase in photon energy in the limit k -» 0 is, from 
equation ( 2 ), 

4K (M) 

A(k -» 0; 0 ) = 1 + - 77 , . (3) 

2 

For 0 « 1 , A -» 1 + 40 , and for © » 1, A -> 160 . These limiting 

forms can also be determined from simpler considerations (Rybicki 
and Lightman 1979). In addition, A(k, 0) = 1 for k= 30, as ex- 
pected since the average photon energy in a Wien distribution is 30 . 

The approximation we employ for Compton scattering is to 
assume that each photon of energy k scatters at a rate defined by 
the exact Compton opacity averaged over a MB distribution of par- 
ticles, K c (k; 0), and is promoted in energy by the factor A(k; 0 ), 
equation (2). Then the source function of equation (1) is given by 

Q c (k,k / ) = cK c (k , )(k'/k)n Y (k') , (4) 

where k and k 7 are related through k = A(k 7 ) • k 7 , and the phase 
space factor k 7 /k is required to preserve total photon number. 
Alternately, introduce the factor p(k) that tells the energy k 7 that 
a photon of energy k had prior to being Compton scattered: 
p(k) • k = k 7 . The static solution to equation ( 1 ) is 


n (k) 
Y 


S (k) K (k 7 ) 

— v + _£ 

cK tot^ k tot^ 


n y (k 7 ) 


(5) 
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where the term in braces 


the LHS of equation (5), n 


with k -» k 7 and k 7 -» k 77 , 


in equation (1) is denoted K^^^k). On 
(k 7 ) is given again by equation (5), but 
so 


n (k) 
Y 


where 


[clC TOT (k)] 


-1 


i = 0 


rf x j (k) 

j=0 


S y [p 1 (k)] , 


X 3 (k) . I P j 00\ 


K ToT [p ;i (k)] 


( 6 ) 


(7) 


The argument p*(k) is an i order operator on k: p (k) = P(k) • k 

= k 7 ; p (k) = p(k 7 ) • p*(k) = p(k 7 ) • k 7 ; and so on. At large optical 

depths, it is only the presence of free-free absorption that shuts off 

higher order scattering and provides convergence, as S (k)~ l/k 

Y 

for bremsstrahlung photons in the limit k -» 0. Hence the importance 

of low frequency absorption and production processes in the limit of 

large optical depths. Equation (6) can be numerically evaluated to 

determine n (k) and z. 

Y 

Because K = K [k, n (k)l , equation (5) must be solved in 
Y YY Y 

an iterative sequence. A steady state solution requires the equality 
between the pair production and pair annihilation rates: 


r (0) 
a 


r 

par -par 


(©) + r (0) + r (0) . 

y -par y -y 


( 8 ) 


Expressions for these quantities in terms of integrals over cross 
sections and photon distribution functions are given in the preceding 
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chapter. Joint satisfaction of equations (6) and (8) yields the photon 
number spectrum and the allowable steady state positron density z. 
Numerical solutions near 0 — ©(1), where simple approximation 
techniques are not possible, will be shown later. We will first con- 
sider features of the results that do not require extensive numerical 
simulation. 

B. Ultrarelativistic Regime: Approximate Treatment of 
Comptonization 

The first systematic study of the physics of relativistic 
thermal plasmas is given in the paper by Bisnovatyi-Kogan, 
Zel'dovich, and Sunyaev (1971). Their principal conclusion is that 
a stable, optically thin, relativistic plasma cannot exist above a cer- 
tain critical temperature © cr . Even in the limit of zero proton 
optical depth in which the photon density goes to zero (Table 1), pair 
production through particle -particle collisions overwhelms pair 
annihilation, and any increase in energy goes into pair production 
rather than thermal energy (Stepney 1983). The existence of the 
critical temperature is due to the decrease with energy of the annihil- 
ation cross section a , equation (III- 2b), and the corresponding 

SL 

increase in the pair production cross sections o , equation (III-27). 

PP 

We can summarize the argument as follows: in the limit of 
zero optical depth, the pair balance relation (8) gives 

„ + n _a a (Y r ) «. <n + +n_)a pp (Y r ) . 


(9) 
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where y is a characteristic electron Lor entz -function at tempera- 
r 

ture 0 . The less important proton-electron pair production cross 

section is neglected (see Figure 9). Solving the quadratic equation (9) 

2 

for n, , and finding that the determinant a - 4 G a vanishes 
+ a a pp 

__ _ 2 
when v = 850 gives © — 17, since v « 3© 

1 r 6 cr r 

The original paper of Bisnovatyi-Kogan et al. gives 0 « 40, 

as does Lightman (1982). Zdziarski (1982) points out that the ther- 
mal averaged pair production rate from particle collisions used in 

these studies is a factor of four too low. He calculates © — 25, 

cr 

so this crude estimate is not far off. Calculating the thermally aver- 
aged pair production rates from particle collisions through equa- 
tions (II- 16) and (111-27), along with the pair annihilation rate (III- 5), 

gives an accurate determination of © .A linear increase in the 

cr 

cross section, as suggested by Svensson (1982), from threshold to 
= 50 for the electron-proton rate, and from threshold to y^ = 100 
for the electron-positron rate, gives the rate coefficients R atl d 
shown in Figure 9. The pair balance equation (8) becomes 

z(z+l)R = (2z + l)R + [(z 2 /2) +(z + l) 2 /2 + z(z+l)]R , (10) 

a ep ee 

dropping all photon- related pair production rates in the limit T 0. 

P 

Equation (10) has no positive real solutions for 0 >0 — 24, using 

cr ^ 

the values shown in Figure 9, thus agreeing with Zdziarski's result. 


Physically, at temperatures above 0 , any incremental 
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Figure 9. Reaction coefficients vs. temperature 0 in relativistic 
plasma. Numerically calculated rate coefficients are 
plotted as a function of temperature, including pair pro- 
duction from particle -particle, photon- par tide, and 
photon-photon interactions, in addition to the annihilation 
rate coefficient. For pair production from bremsstrahlung 
photons, only the rates associated with the electron- 
proton bremsstrahlung spectrum are shown. Symbols for 
rate coefficients are related to rates explained in text. 
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increase in the pair density will be accompanied by an increase in 
the rate of pair production through particle -particle collision that is 
always greater than the corresponding increase in the rate of pair 
annihilation. Therefore no pair balance is possible at these tem- 
peratures. If the optical depth of the system is non- negligible, 
internally produced photons will also contribute to pair production. 
One should think that the effect of this would be to lower © and 


cr 

give a critical temperature ® cr J T p)- This is indeed the case, as 
investigated by Lightman (1982) at extreme relativistic (ER) temper- 
atures. In fact, © = © (r , n ), from preceding considerations. 

cr cr p p 

cr 

In what follows, we calculate © , or equivalently T (© , n ), in 

cr p p 

the limit © » 1, in addition to the steady state positron density z, 
on the basis of approximations to equation (6). 

At temperatures greater than about 2 MeV, that is, for 0^4, 
asymptotic expressions for the rate functions in equation (8) can be 
evaluated analytically. For those rates depending on the photon dis- 
tribution, it is necessary to determine an approximate form for the 
photon number spectrum. This we do by assuming that the photon 
distribution is composed of two components in the pair production 
region: the internally produced bremsstrahlung photons and the 
Compton scattered photons. We further assume that those photons 
Compton scattering into the regime of photon energies k <M@) 


obtain the Wien form 



100 


w 


w 

dn (k) = 

Y 


n 


2 ©' 


exp(-|ik)k dk , 


(ID 


W 3 

or f (k) = exp(-|_tk), with A = n /8tt© from equation (III- 26). The 

Y Y Y 

W 

volume density of photons in the Wien distribution is then just n 

The Wien spectrum is the most probable distribution that a 
gas of photons will obtain after Compton scattering with a thermal 
distribution of electrons (Synge 1957). It corresponds to the Bose 
distribution with negative chemical potential, reflecting conservation 
of photon number. One can formally show on the basis of the formal- 
ism developed in Section II. E that the time scale for energy exchange 
through photon-electron (Compton) scattering is always much shorter 
than the time scale for energy loss through bremsstrahlung due to 
the different orders of the two processes; thus, the model is con- 
sistent. However, it would be necessary to solve a time -dependent 
Boltzmann equation using the exact Compton scattering function, 
equation (11-52), to characterize the establishment of a Wien distribu- 
tion after one, or a succession, of scatterings of a photon in the 
region k~ ©(©). The assumption adopted, namely that this occurs 
after a single scattering into this regime of photon energy, remains 
uncertain. 

In the high temperature limit the spectrum of annihilation 
photons gives only a small contribution to the total photon number 
spectrum; for this reason it is neglected. This is because 
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annihilation photons are produced conjointly with bremsstrahlung 
photons, which always dominate at ER temperatures . For example, 
at temperature 0 = 3, the value of the e -e bremsstrahlung spec- 
trum is a factor of 3 greater than the e^-e annihilation spectrum 
at the peak of the annihilation spectrum, as can be seen from a com- 
parison of equations (II-8) and (11-10) and the values shown in Fig- 
ure 3b. Adding in the accompanying e p, e p, and e -e brems- 
strahlung- produced photons, this represents a 15% contribution at 

most, and decreases on either side of the annihilation peak (note the 
2 

k dependence of the Gaunt factor for k» 0 , © » 1 in Figure 3). 
For temperatures © greater than 3, the annihilation spectrum is 
relatively even weaker. 

A corollary to this observation is that no annihilation line 

should be visible from a stable ER thermal plasma (Lightman 1982). 

Indeed, because of the decline of the pair density at subrelativistic 

temperatures, except when the optical depth is so great than the 

annihilation line then suffers severe Compton broadening, Svens son 

(1984) concludes that an annihilation feature cannot be observed from 

a stable relativistic plasma in pair balance. 

In order to calculate the number of photons being scattered 

W 

into the Wien regime per second, S (0), an asymptotic form for the 
electron- proton brems strahlung Gaunt factor is useful. At ©» 1 


and k « 0 , one obtains 
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g(k; © ) 



( 12 ) 


using the approximate form for the electron-proton bremsstrahlung 
cross section, 


da 




dk 




1_ 

2 


(13) 


valid for soft photon production at relativistic (y » 1) electron 
energies. At ER temperatures, the various thermal bremsstrah- 
lung production spectra have identical form. Incorporating appro- 
priate multiplicative factors, the total bremsstrahlung production 
spectrum becomes 


S b (k;0 ) = (8z 2 + IOz + 3) 
V 



exp(- |ik) S/n 



(14) 


for k « © and 0 » 1, where T| = 2 exp(5/2 - 2y ) . 

hi hj 

W 

The value of S (©) is determined from the relation 


W 

S (0) = 


dk p(k; y) S x (k; © ) 


(15) 


where p(k; y) is the probability of a photon with energy k being 
scattered into the Wien regime, and y represents those parame- 
ters (0 , n , . . . ) upon which p may depend. In the Compton scat- 
P 

tering approximation we use, 
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p(k; y) = X(k) a(k) , (16) 

where a(k) = 2n(3® /k ) /bn A(k; ©) characterizes the average number 
of scatterings required for a photon of energy k to reach the Wien 
(k fa 30) regime, and X represents an averaged probability that the 
photon will be Compton scattered rather than absorbed or lost from 
the system. Its meaning will become clear as specific regimes of 
photon energy are treated. 

We require a photon to scatter at least one time to enter the 
Wien regime. It is therefore appropriate to divide the bremsstrah- 
lung generated photons into two regimes depending on whether the 
energy k of the photon is less than or greater than 30/A, where A 
denotes the asymptotic value A (k -» 0; © ), equation (3). This is be- 
cause A is sufficiently large that A (k = 30 /A; 0) » A for 0 » 1. 
In the former regime, a photon has a probability of scattering into 
the Wien regime given by the multiple scattering form of p, equa- 
tion (16). In the latter regime p becomes X(k), the probability of 
a single Compton scattering at photon energy k. In addition, photons 
below a_ characteristic, energy k^ _ defined through the relation 
K^(k^) = K^k^) more likely be absorbed through the inverse 

bremsstrahlung process than lost from the system because of the 
rapid increase in the free -free absorption coefficient as k^ -» 0 
(equation (III- 32)). (In any event, the value of k^ is not critical at 
low to moderate optical depths because these low energy photons 
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will only rarely saturate to the Wien regime. ) Equation (15) be- 
comes 

/ 30/A 

dk X ai ; (k) S (k; 0) 

Y 

k b 

00 

+ / dk X(k) S^(k; © ) . (17) 

*3© /A 


An examination of the absorption coefficients in Figures 8 

and 10 shows that the only significant sources of photon opacity at 

photon energies k^ < k < 30/A are loss from the system and 

Compton scattering. The Compton scattering opacity for these 

energies takes place almost entirely in the Thomson limit. This is 

because the separation between the Thomson and Klein-Nishina (KN) 

regimes is at photon energy k determined by the condition y 

s r 

(k « 1), where “ is again the Lorentz. factor typical of thermal 
s r 

electrons. For 0» 1, we have k » (3®) ^ . This may be com- 

s 

pared to the value k = (3© /A) (50) * separating the two scatter- 

ing regimes in the ER limit. 

Therefore, when k^ < k < 30/A, we can calculate X(k) by 
assuming that the scattering is described in the Thomson approxima- 
tion to the Compton cross section. Hence X(k) — K . (K + K ) 

— [l + (t ° + T ° 2 /3) = Xq in sf (0), where T ° = T c (k -» 0) = 

(2z + 1)t . The limits of X are X -> T° for « 1, and 

p U U O L/ 
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2 ^ 

X n -» 1 - (3t ° )” for T° » 1. Integrating S , using equation 
0 C/ i 

(14), one obtains 


W 

! i» = I b 


30 /A 
/ ■ 


, -1 Jn(3®/k)/&i A _2 

dk k X v ^/i(ri 0 /k) 

0 -k 




(18) 


E • Q • i Q 
b 


' U 1 / T1 E®\ % I 

X o\ hX 0- 1 >-V X 0- 1 rm (x 0 - X 0>( • 


A „ J 


where Q = Z/nA/S/n X , u,-, = S/n(?>®/)<L )l%/n A, and E, = 

0 r> r> D 

2 2 

(8z + lOz + 3). (2n aca_/TT). In the case of unsaturated Compton- 

P T 

ization of the low energy photons k > k^ » that is, when X’ « 1 , 


Sj ( 9 ) « S b - o • X 0 [Q(1 -fi«x 0 ) - x Q a« (n E ®/ 3 ] . (19) 

w 

We now turn to the second integral of equation (17), (©), 

representing those photons produced with k> 3© /A which scatter 
one time before entering the Wien distribution. It is necessary to 
incorporate the effects of the KN decline in the Compton cross sec- 
tion on the Compton opacity in order to calculate this rate. At ER 
temperatures, the Compton opacity is given by Gould (1982a)_or from 
equation (III- 24b), namely 

K c (k;®) = (2z+l) ( 3/ 1 6 k0) • S/n(r\ ^ k© ) , (20) 

where T)^, = 4 exp(l/2 - y ) . This expression is valid for 
k> 1/©, so we need an interpolation formula between the Thomson 
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and KN regimes to account for photon energies in the range 3®/A < 
k < l/@ . Take k = 30 0/A for the upper photon energy of the inter- 
mediate regime, which is well within the KN regime of the Compton 
opacity. One finds that K (k = 300/A) — 0. 19 (2z+l) n for 

C p JL 

0 » 1. An accurate (10-20%) interpolation formula at 30/A ^ k < 
30®/A is given by the formula 

K (k; 0) = (2z+l) n CT_ (1 - 0. 43 k©) , (21) 

p l 


which reduces to the known values at k = 30© /A and for k -» 0. 

Inserting equations (20) and (21) into the second integral of 
equation (17) gives 


s S b (2a+1)T P 


.30 0 /A 


(1+t“/3) 


dk k _1 (l 


0. 43 k©) 071 


3 0 /A 



+ 



071 (r| k©) Bn 



( 22 ) 


in the approximation T^. < 1. The exponential factor from equa- 
tion (14) is taken equal to unity in the first integral above, a good 
approximation considering the photon energies of the intermediate 
regime, and is used to impose an artificial cutoff at k = 3® in the 
second integral. This is a crude approximation, not to mention the 
fact that the Gaunt factor (12) is no longer valid at these photon 
energies. But the contribution of the second integral in equation 
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(22) is at any rate small compared to both and the first integral 

in equation (22) at ER temperatures because of the secular decrease 
in the Compton opacity (20). 

Evaluating the integrals in equation (22), retaining only the 
dominant terms from the second integral, gives 

5^ = Z (2z+l ) t [(1 + t°/3) • 4.7 &»( 2. 50) + 0. 1 fln(ri _0)«n(ri „0 2 ) ] , (23) 

b D P Vj V^/ -E_I 

2 

using A = 16© in the limit ©» 1. This equation, in conjunction 
with equation (18) or (19), provides an approximate analytic treatment 
for the rate at which bremsstrahlung photons are Compton scattered 
into the Wien component of the photon number spectrum per unit 
volume. Though equation (23) is only valid under the assumption 
Tp(2z+1) < 1, numerical comparisons of equations (18) and (23) show 
that when t° > 0.1, that is, for low to moderate optical depths, 
equation (23) gives a small contribution to equation (17). This simply 
reflects the cascade of low energy photons into the Wien distribution 
as the optical depth for Compton scattering begins to be significant. 
Since equation (18) does not suffer from the low optical depth approxi- 
mation, this analysis permits a treatment of Comptonization in an 
ultra -relativistic plasma at all optical depths. For an alternate treat- 
ment of Comptonization in a thermal plasma at ER temperatures, the 
reader is referred to the paper of Lightman and Band (1981). 
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Figure 10. Photon-particle pair production opacity. The opacity 
of a photon with energy k due to photon- particle pair pro- 
duction from a MB distribution of electrons or protons is 
shown at various temperatures. For photon- photon pair 
production at these temperatures, the opacity is given by 

n ct (k), where n is the proton number density and 
P YP P 

CT (k) is the cross section (II-Z9). 

YP 
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C. Untrarelativistic Regime: Asymptotic Forms for the 
Rate Processes 

Having described a method for calculating the rate at which 
photons populate the Wien component, it is now necessary to give ex- 
pressions for the rates associated with pair production and annihila- 
tion from the various photon and particle processes. At ER tempera- 
tures, asymptotic forms can be derived so that the pair balance equa- 
tion (8) can be solved directly. This ultimately gives the steady state 

positron density z(t n ) and the critical proton optical depths 

P P 
c r 

T (n , 0), when equations (8) and (15) are jointly satisfied. 

P P 

The asymptotic expressions given below for rates of pair pro- 
duction through different processes are accurate to roughly 10-20%, 
with the accuracy typically increasing with temperature. This is 
adequate since, as mentioned earlier, the basic assumption that the 
particles obtain a thermal form is not correct at these temperatures 
(Section III. F), though an effective temperature could be assigned on 
the basis of the average particle energy. The annihilation rate r 

a 

at ultrarelativistic temperatures has already been derived. It is pro- 
portional to the reciprocal of equation (III- 42): ' “ 


2 2 2 
T a - z(z+l) n (c TTr^ ) ® /(2® ) 


(24) 


valid for 0 »1. The associated rate coefficient is shown in Fig- 


ure 9. 
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For the production of pairs in e -e collisions, the reaction 


rate 


2 2 3 2 

r, = z(z + l) n 28 a o T c to 0/9n 
+ - p T 


(25) 


is used (Zdziarski 1982). Since pair production from e -e and 

+ + 

e -e collisions have similar form, the total pair production rate 

r for the three mechanisms is equal to r , , but with z(z+l) 

ee + - 

2 

changed to (2z + 2z + l/2). An adequate representation for the pair 
production rate from collisions of protons with electrons and posi- 
trons is, fitting the numerical results shown in Figure 9, 


r = (2z+l)n^ c 2 ti^©/16. (26) 

ep p T 

The rate of pair production through photon- particle and 

* 

photon-photon interactions depends on the steady state photon spec- 
tral density. For the internally produced bremsstrahlung photons, 
the spectral density is given by 

n b (k) = T t S b (k) (27) 

Y P c Y 

where t„ = (c n ct„) ^ and S b (k) is the source function of brems - 
v> pi. y 

strahlung photons. ' This relation assumes that the magnitudes of the 

b i b 

absorption coefficients K _ , K , and K are small com- 

C/ Y -par Y ■ Y 

pared to K at photon energies important for pair production, that 
is, for k> 1/0 (cf. equations (1) and (5)). 

In other words, the use of equation (27) to calculate rates 
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from a bremsstrahlung-type spectrum of photons does not require the 
system to be optically thin except at photon energies k > 1/0 . Other- 
wise there is a frequency-dependent distortion, and the use of equa- 
tion (27) will not provide valid estimates for the pair production rate 
processes. In terms of the Compton opacity, this requirement trans- 
lates into = K (k) • L/3 at k = 1/0 , where the Compton opac- 

ity is largest. At ER temperatures, this value of k is in the KN 

regime and one gets the condition T° - (2z+1)t <12, using equation 

L. p 

(20) to evaluate T^(k). 

Absorption of energetic photons through pair producing mech- 
anisms can also distort the spectrum (27). From equation (23) and 

the high energy form of the cross section (III- 28), o (k) -» 

ye 

(7 a /6 tt)ct t (0'/i 2k - 109/42), one obtains 


K (k) = (2z + l)(an ct t ) 
y e pi 





&»( 25 ) - ^ 


+ e U5/2 E 1 (u6/2)[ , (28) 


where 6 = max[6 + , w/2] and 6 + — 4/k. For 4«k«80, equa- 
tion (28) reduces to 

1C (k) = (2z+l)(7 /6 tt)( an o ) !/n (4 r\ k© ) , (29) 

ye p r p 


with r\ =exp(l^y r ,- 109 /42). 
P ^ 

at photon energies where K (k) 


Requiring K < n (J„/t 

ye p T p 

is largest in the pair producing 
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regime, namely for k~ 3© , one obtains the requirement 

(2z+l) t < 6TTa ^/7 bn(4r\ k©), a much weaker condition than that 
P P 

imposed by the Compton opacity. The contribution to absorption due 
to pair production in photon- proton collisions is even smaller than 
K (k), though the asymptotic forms of the cross sections (11-28) and 
(111-29) are the same in the limit k -> 00 . The argument resembles 
the one presented in the Appendix contrasting reaction rates for a 
nonrelativistic and ER system of particles. 

To determine whether the bremsstrahlung spectrum is dis- 
torted from photon- photon absorption of one bremsstrahlung photon 

b 

by another, one must calculate the absorption coefficient K of a 

YY 

photon traversing a photon gas having a bremsstrahlung -type dis- 
tribution. Using the spectrum (14) for the bremsstrahlung photons 
with equation (III- 25) gives 


K b (k) = T (8z 2 + 10z + 
YY P 


2 r -3 /^E® k \ _ 

3)(2n^ar e ) I duu bn I - I cp (u) . (30) 


The values of the definite integrals defined by equation (30) are 1. 55 
and 3. 40 from numerical calculation, the latter number referring to 
the integral containing the bn u term. From equation (30), one ob- 
tains the condition 


(8z 2 + lOz + 3) t 2 < a ^ /l. 2 bn(r\ k® 2 ) , 
p W 


(31) 
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where ri • = 0.86, requiring that K (k) < n a_/T . This esti- 
W YY p T p 

mate is accurate because pair production in photon- photon collisions 
occurs chiefly for photons k and k* related by k • k ; ~ ©(1) due to 
the threshold and high energy behavior of the photon- photon pair pro- 
duction cross section. When k~ ©(*«)), the bremsstrahlung photons 
k / ~ ©(9 1) primarily responsible for pair production are well de- 
scribed by the Gaunt factor (12) at ER temperatures. 

From the preceding considerations, we see that the severest 
restriction to the use of the spectrum (27) for calculating rates comes 
from the Compton opacity condition, T° < 12. Therefore, at optical 
depths as great as one or two, the calculated pair production rates 
will be accurate. In fact, when optical depths approach this value, 
the greater portion of pair production proceeds through the photons 
Comptonized to the Wien form, so that the rates calculated from 
equation (27) are least accurate when the net contribution of the 
bremsstrahlung spectrum (27) to pair production is smallest. 

Numerical techniques were used to evaluate pair production 
rates from a bremsstrahlung spectrum containing the exact B-H 
formula as a check to analytical estimates. The formula derived by 
Lightman (1982) for the pair production rate from photon-electron 
collisions accurately describes the numerical results for this process. 
The rate is 

r y e = T p(^ z+ l)(8 z 10z+ 3) (~-) ( a^n^co^) 2/r?© Urn. 4® (1 - 0 *) , (32) 
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adding the term (1-0 *) to the equation of Zdziarski (1982). The 
expression 

b 2 2 2 2 3 

r = T (8z + lOz + 3)(0. 014)(a n co_) Bn 0 Q/n 30 (33) 

YP P P T 

is valid for the pair production rate due to brems strahlung photons 
interacting with protons at ER temperatures. The pair production 
rate due to self-interacting brems strahlung photons is described by 
the formula 


b22 2 22 3 2 

r = T (8z + lOz + 3) (0. 91)( a n ca T )to 0(1 + 2/0 ) . (34) 

YY P P T 

The numerical calculations for the associated rate coefficients are 
shown in Figure 9. 

The values calculated for r are in poor agreement with 

YY 

Zdziarski's numerical calculation. Moreover, the rates given by 
Svensson (1982) for these processes are seriously in error, so the 
results of the calculations of this paper are suspect. 

Turning now to the Wien component of the photon spectrum, 
considerations based on equation (1) give 


W 

n 

Y 


[ c(K + K + K )] _1 S W 

L y-par yy y 


(35) 


for the number density of photons in the Wien distribution. The pro- 
duction rate of Wien photons was calculated in the preceding section, 
and the term in brackets represents an average residence time 
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— W 

T (®) for a photon in the Wien distribution ( K , of course, is 
absent, since Comptonization merely sustains a Wien photon in the 
distribution). Letting 


T W (@) = [T“ 1 (@) + T' 1 (0) +T‘ 1 (0)]‘ 1 f 

L y - par YY 


(36) 


it is a simple matter to calculate the individual average timescales of 
equation (36) by considering averages over the Wien spectrum (11). 
Recalling the form of K (k), we have 


oc 

T (0) = (2© 3 ) 1 / dk k 2 (c K^) 1 exp(-ptk) 

•'A 


C 

'l 


= t c T (2© 3 )" 1 I dk k 2 (l +t c / 3) exp(-Kk). (37) 


For © » 1, an insignificant fraction ~(6®^) * of the Wien photons 
are scattering outside the KN regime k> 0 . Thus, it is an excel- 

lent approximation to use equation (20) in equation (37), giving 


T (0) = t T 11 + 

p 


o 2 

T c 0n (r) L 0 ) 

320 2 


(38) 


where r) =4 exp(3/2 - Zy ). 

Li . ilj 

There are four other average residence times in equation (34) 
associated with photon-electron, photon- proton, and photon- photon 
absorption of the Wien photons, the last depending on whether the 
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opacity results from photons in the Wien or bremsstrahlung compo- 
nent of the photon spectrum. These timescales may be estimated as 
proportional to the inverse of the reaction rates, according to the 

WWW w 

relation T = n /r (0), with analogous expressions for T , 
yeyye yp 

W Wb 

T , and T , employing obvious notation. 

y Y y Y 

The particle production rates from photon-particle collisions 
can be derived from the asymptotic forms of the cross sections 
(III- 28 ) and (III-29) at high energies. One obtains for the pair pro- 
duction rate from ER thermal electrons and positrons interacting 
with photons in the Wien distribution, using equations (III- 30) and 

(29), 

r^ e (@) = (2z+l) aca T n p n^ (y^y) ^(r^© 2 ) , (39) 

where T) = 4 exp(5/2 - 2y^, - 109/42). The pair production rate 
from Wien photons interacting with protons can also be determined 
from equation (III-30), though now with (i >> 1. The result is 

W ^ W / 7 

r ( 0 ) = aca T n n 77 — 
yp T p y \16 tt 

The term in © * was added in order to more accurately fit the 
numerically calculated rate near 0 =: 3 . 

The pair production rate from a self-interacting Wien dis- 
tribution of photons can be derived from equations (III- 31) and the 
asymptotic form of cp (k * k') given by Gould and Schreder (1967). 


) (0/10 - 0. 98 + 1. 2/0) . (40) 


One obtains 
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W 

r 

YY 


( 0 ) 


2 W 2 
cTTr n 
e y 


[0/1 © + (Bn 4 - 


2y e )]/2© 2 


(41) 


at ER temperatures. The pair production rate from the Wien com- 
ponent interacting with the bremsstrahlung component of photons is 
calculated from equations (III- 30) and (30). One finds that at ER 
temperatures, 


Wb 2 2 W 3 

r ( 0 ) = T (8z +10z + 3)3. lacr nn 2n( T| ® ) , 

YY p e p y Wb 


(42) 


where r l-^ r ^ = T l-y^ exp(3/2 - Yg )• This expression is accurate 
to 5% for temperatures © > 3, as verified by the numerical calcula- 
tions plotted in Figure 9. 

The preceding rates were used to model a stable ER thermal 

plasma. The results are shown in Figure 11. In the Figure, the 

positron density is plotted as a function of temperature for various 

values of the proton optical depth , at proton density n = 

10 -3 

10 cm . Either zero, one, or two roots exist for each value of © 
and Tp , because the pair balance equation (8) is a polynomial equa- 
tion in z. The annihilation term, proportional to z(z+l), vanishes 
at zero, whereas the pair production term is a higher order poly- 
nomial with positive coefficients multiplying the powers of z, so the 
maximum number of roots is two. 

The physical conditions corresponding to the absence of a 


stable solution to the pair balance equation were discussed above. 
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Figure 11. Equilibrium, pair density and maximum proton optical 

depth at extreme relativistic temperatures. The steady 

state pair density z of a nonmagnetic relativistic thermal 

plasma is shown in Figure 11a as a function of temperature 

© for © » 1 at different values of the proton optical 

10 - 3 

depth T . The proton number density n^ = 10 cm 

The maximum t of a steady thermal relativistic plasma 
P 

is shown in Figure lib as a function of critical temperature 
© cr . For proton optical depths greater than the values indi- 
cated, a stable system must be in a condition of complete 
thermodynamic equilibrium. The dependence on proton 

number density at these temperatures is shown for 

-3 25 -3 

n =1.0 cm and n =10 cm 

p p 
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The temperature at which a single root exists represents the maxi- 
mum critical temperature for a given proton optical depth. Values 

of T (0 ) at ER temperatures are shown in Figure lib. The 

p cr 

maximum critical temperature is approximately equal to 24 in the 
limit of zero optical depth, as calculated earlier. A weak depend- 
ence of the critical temperature on n^ is observed at densities of 
astrophysical interests, which becomes more pronounced as in- 

creases, reflecting the increasing importance of Comptonization and 
the value of the bremsstrahlung absorption cutoff. The model breaks 
down at temperatures 0 ~ 2-4 where the asymptotic forms for the 
rate factors begin to fail. Discussion of the double roots below the 
critical temperature at a given optical depth will follow a treatment 
of the system at subrelativistic temperatures, where the same be- 
havior is observed. 

D. Subrelativistic Regime: Large Optical Depths 

For temperatures 0 « 1, the steady state positron density 
decreases exponentially with temperature at low optical depth due to 
the threshold behavior of the cross section and the Boltzmann factor 
exp(-E/0 ) governing the distribution of particles or photons with 
energy E. This behavior is suggested by the variation of the rate 
coefficients in Figure 9 at low temperatures. One sees that a sig- 
nificant positron density will be produced only if there is a large 
Wien component of the photon spectrum, that is, when the system 
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has large optical depth and Comptonization is important. The large 
W 

value of R ^ compared to the other pair production coefficients 
obtains because the average photon energy in the Wien distribution 
is relatively large. 

The phase diagram corresponding to Figure 11a can be com- 
pleted at subrelativistic temperatures and large optical depths by 
applying the method of the last two sections to this regime. Compton- 
ization will be treated as before, though now with the appropriate 
production spectra. The most important photon sources are annihila- 
tion, bremsstrahlung, and double Compton production. The double 
Compton mechanism if of the same order as bremsstrahlung in the 
fine structure constant, so is significant when the density of photons 
approaches the particle density. 

We employ the nonr elativistic Born approximation for the 
bremsstrahlung production cross section in the soft photon limit at 
subrelativistic temperatures, giving 


b -i /V 

S (k; 0) = E k fa -f— 
y b \ k 


(43) 


where S = 2[z(z + l ) + 2^ 2 (1 + 2z)] n ^ a ca„( tt®) and r\ = 

b P T g 

4 exp(-Y-g). The term proportional to z(z + l) accounts for the con- 
tribution from e -e bremsstrahlung, and follows from the consid- 


erations of Joseph and Rohrlich (1958; cf. Appendix). The other 


+ + 


term is for e-p bremsstrahlung. The e -e and e -e thermal 
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bremsstrahlung spectra are much smaller at these temperatures 
since there is no dipole moment in the scattering system of two 
identical particles, so is neglected. In the same soft photon low 
temperature regime, the double Compton spectrum 

(k; © ) = (2z + l) n n^ ac r ^ ('7~) ® k (44) 

Y p y e \ 9 / 

is used (Svensson 1984). 

Since © « 1, all Compton scattering takes place in the 
Thomson limit. Solving the integral (15) for the rate of population 
of the Wien regime by Comptonization, noting that for ^ 1 most 

photons originate in the soft photon portion of the production spectra, 
one obtains 


TOT 


r a<9) + E dc +*.H 




? r \ 
+ s b n x o 


_ X () b <u b *0 ' M ‘ X o° (u 0 x 0 ‘ 1 _ 


The terms X , Q, and u, are defined as before, and u. = 3 /^/i A , 

u b 0 

though now A -> 1 +40 in the definitions of u^ and u^ . The value 
of k^ is found by identifying K ^(k) and K^,(k), since K « at 
large optical depths. The analytic result (30) was evaluated by im- 
posing a cutoff at k = 0 which, from our previous discussion, should 
not be a bad approximation in an optically thick system. 

The first term in equation (45) is due to the Comptonized 
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annihilation photons, the factor of two accounting for the two photons 
produced per annihilation event. At temperatures © ^ 0.4, the 
annihilation rate is given approximately by 


r (©) = z (z + 1 )n^ R (®) = z (z+l)n ^ cTTr ^ (1 - 0. 9® ) . (46 ) 

a pa p e 


The probability that an annihilation photon will be Comptonized rather 
than ultimately absorbed or lost from the system is given by the ratio 
of the absorption coefficients 


1 + 


( K 


K 


+ K a + K Wa 
L YY YY 
K „ 


-1 


TOT 


( 47 ) 


neglecting the opacity that an annihilation photon might suffer due to 
pair production from a proton, positron, or electron, which are expo- 
nentially small at these temperatures. The absorption coefficients 
Wa a 

K and represent the opacity of an annihilation photon to pair 

production from interactions with other photons in the Wien or annihil- 
ation component, respectively. At these low temperatures most 

annihilation photons have energies near k=l, so we can approximate 

Wa 

the coefficients by the frequency- independent relations cK = 

YY 

2 , cL 3, cL cL 

r (®)/n and cK = r (®)/n . where n is the number density 

YY Y YY YY Y Y 

of annihilation photons. Numerical expressions for the rate coeffi- 

Wa a 

cients r (©) and r (©) to photon-photon pair production are 
YY YY 


given below. 



124 


W 

The number density of Wien photons n is computed as in 


Y 


— W 


the earlier section, by introducing an average residence time T (©) 
during which a Wien photon survives before being lost from the dis- 
tribution. It is given by 

-1 


— W 

T ( 0 ) = 


x W _ . W x Wa, m , W 
cK + r (0)/n +r (®)/n 
L YY Y YY - Y 


(48) 


W 


retaining only the important absorption processes. Thus n 

Y 

— w w 

T (0)S (0), with the expression for Wien photon /annihilation photon 

Wa 

reaction rate r given below. 

YY 

The reaction coefficient for the production of pairs due to a 
self- interacting Wien distribution was numerically calculated and is 
shown in Figure 9. Weaver (1976) gives a numerical'fit to the low 
temperature rate. Noting the errors in his equation (33), it is 


2 2 - 2/0 

cTT r e . yy2 

r = — (1. 042) [1 + 0. 7280 ] ' n 


YY 


8 ©' 


Y 


(49) 


good to 5% for 0 < 0.4. A numerical fit to the rate coefficient 
Wa 

R (0) describing pair production from annihilation photons inter- 
YY 

acting with Wien photons is given by 


R Wa (0) = 0.80 cr 2 ©” 1 * 9 e' 1/@ 
yy e 


(50) 


The associated reaction coefficient is given by the expression 


r Wa (0) = z(z + l) n 2 n W R Wa (0) t-, 1 T a (0) 
YY P Y YY c 


( 51 ) 
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a 

employing the value t^, for the average lifetime T (0) of an annihil- 
ation photon in the calculation of equation (51). Likewise character- 
izing the average spectral density of annihilation photons through the 
use of the residence time t^ , one obtains 

r a (0) = z^(z + l) n^ [2.06 cr^0 exp(-40)l t_^ [T a (0)] (52) 

YY p e C 

for the rate at which pairs are produced by self- interacting annihila- 
tion photons, from a fit to the numerical results. The rate coeffi- 

3 . 

cient R is the term in brackets in equation (52). Equations (50) 
and (52) are accurate to 10% for 0. 075 ^ ® < 0.4. 

““cl 

The mean lifetime of an annihilation photon is given by T = 
(c^TOt) 1 * t ^ ie definition (47) of and the rates (46), 

(51), and (52), one finds that 

-1 


(53) 

Given n , T , and z, equation (48) and (53) are iterated until the two 
P P 

are jointly satisfied. Then z is incremented until the pair balance 
equation 

a W a Wa 

r (0) = r (®) + r (0) + r a (0) (54) 

yy yy yy 


-1 -a 


-1 


t r T (0) = \ (2z+ 1 ) + [ T (1+T / 3 )] + 


r , . 2 a W W a , 

Lz(z + l)n R + n R i 

p yy y yy 

2R 


is additionally satisfied. 
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The numerical simulation was performed throughout the tem- 
peature regime 0 < 0. 4, not only where the expressions for the rates 
associated with annihilation photons were accurate. This was be- 
cause these rates proved small at the temperatures and optical depths 
used, relative to the pair production rate (49) associated with self- 
interacting Wien photons, which is valid down to the lowest tempera- 
tures. 

The results are plotted in Figure 12. for a background proton 
10 - 3 

density n = 10 cm , along with the ER temperature results for 
comparison. Proton optical depths T =0.5, 1.0, 3.0, and 5.0 
were used. Although the approximations (such as the calculation of 
k^) applied to the optically thick regime, it appears that one could 
smoothly extend the curve z (0 , T = 0.5) to the expected values at 
ER temperatures. Additional pair production due to internally pro- 
duced bremsstrahlung photons would serve to increase the value of z 

in the low z branch of the curve for T =0.5, improving the agree- 

P 

ment. Thus we can be confident that the results for the model at ER 
temperatures and at subrelativistic temperatures, large optical 
depths are consistent. 

E. Discus sion 

Above a critical temperature, a steady optically thin non- 
magnetic plasma cannot exist. The efficiency of particle producing 
mechanisms will ensure that excess energy is transformed into 
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Figure 


12. Equilibrium pair density at subrelativistic tempera- 
tures. The steady state pair density z of a nonmagnetic 
relativistic thermal plasma is shown as a function of ® 

for ® « 1 at different values of the proton optical depth 

10 - 3 

T . The proton number density n =10 cm . The 
P P 

results at extreme relativistic temperatures, from 


Figure 11a, are shown for comparison. 
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particles and photons until conditions of complete thermodynamic 
equilibrium obtain. To produce a system in thermodynamic equilib- 
rium at relativistic temperatures would require a tremendous energy 

4 

source, as the energy density goes as © . 

This conclusion does not invalidate the existence of relativis- 
tic nucleonic gases and the calculations of Section III. E, provided 
that the system has a transient lifetime. For a consistent character- 
ization, the system lifetime would have to be bracketed between the 
timescale for thermalization and the timescale for particle produc- 
tion. The most obvious candidate for such a configuration would be 
accretion into a massive compact object or black hole. Zdziarski 
(1982) has considered these timescale in the context of spherical 
accretion, and a great deal of effort is being devoted to the task of 
calculating spectra resulting from disk and spherical accretion 
(Eilek 1980; Eilek and Kafatos 1983; Meszaros and Ostriker 1983; 
Meszaros 1983). 

The existence of two roots for the positron density below the 
critical temperature at a given set of values of the system parameters 
follows from the assumptions adopted for the thermal plasma. The 
numerical results show that except near the critical temperature and 
when optical depths T > 0. 1, the lower branch is dominated by pair 
production from particle-particle collisions, and the upper branch is 
dominated by pair production through photon processes. When 
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T > 0. 1, it is the Comptonized Wien photons that make the largest 
contribution to pair production in both branches. 

Both solutions are valid for a uniform system. Ldghtman 
(1982) interprets the upper branch as being physically unstable to 
density and temperature perturbations. He arrives at this conclusion 
because the upper branch has a negative specific heat, that is, a dif- 
ferential increase in the energy of the system will lead to a reduction 
in temperature, because most of the added energy will be trans- 
formed into particle mass energy which in turn must be thermalized. 
The physical result will be clumping in the system. 

Comparing the results of the positron density phase diagrams. 
Figures 11 and 12, with the work of other investigators, one finds 
acceptable agreement with only one other study, that of Svens son 
(1984). At the highest temperatures the phase portrait differs from 
that of Lightman (1982) due to his error in the pair production rate 
from particle-particle collisions. At temperatures near 0=3 and 
optical depths at which this mechanism is no longer important, gen- 
eral agreement is found. Disagreement with the paper of Svens son 
(1982) was traced to numerical errors in his calculations of pair pro- 
duction rates arising from internally produced bremsstrahlung pho- 
tons. Evidently these errors are corrected in his 1984 paper, as a 
different phase portrait is presented. Even though different methods 
were used, the numerical values of z(0, T , n ) between this work 

9 ’ rx' -rv 
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and Svensson's latest results are in agreement (20-40% difference at 
worst). 

A complete numerical simulation of the system was also under- 
taken with the Compton scattering process performed by means of 
numerical calculation of equations (6) and (7). It was found necessary 
to broaden the Comptonized spectra over one decade in energy con- 
sistent with conservation of energy and photon number in order to re- 
duce numerical artifacts originating from the discrete nature of the 
Compton scattering approximation. At any rate, this broadening 
probably more closely simulates actual Comptonization by relativistic 
thermal electrons, judging from the Monte Carlo calculations of 
Pozdnyakov, Sobol 1 , and Syunyaev (1978). 

Calculation of the emergent spectra from relativistic thermal 

plasma at temperatures 0 = 1 and 0 = 3 is shown in Figure 13 for 

arbitrary values of the proton optical depth. The system is in a 

10 3 

condition of pair equilibrium with n^ = 10 cm . The lower curves 
of Figure 13 labeled "1" represent the production spectra due to 
bremsstrahlung and annihilation processes. The upper curves labeled 
"2" represent the production spectra after suffering effects of absorp- 
tion, and the curves labeled "3" are the emergent spectra after account- 
ing for Comptonization of the produced photons, calculated in a self- 
consistent manner. One observes moderate Comptonization in Fig- 
ure 13a, and the formation of a Wien peak in Figure 13b. In Figure 
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13c, a distorted Wien hump is seen in the emergent spectrum. This 
is because the temperature used is very close to the critical tempera- 
ture at this optical depth, where photon-photon absorption becomes 
important for photons with energies near k = 1. The features of the 
emergent spectrum at low photon energies in Figure 13c are not to be 
considered as artifacts of the Comptonization method. Comparing 
the Monte Carlo results of Pozdnyakov et al. just cited, one finds 
that low energy spectral variations can propagate over many decades 
of photon energy, even after multiple Compton scatterings. 

A recent preprint (Zdziarski 1984) gives spectra from a full 
Monte Carlo simulation of a nonmagnetic relativistic thermal plasma, 
even treating changes in photon density over the spherical region. 
Though the parameters used are not the same, one can compare the 
features of the spectra in Figure 13 with those of this paper. The 
agreement is good. One must note, however, that Zdziarski arti- 
ficially imposes a pair density in his model, so his system is not in 
pair balance. 

The fully numerical treatment gives general agreement with 
results for the equilibrium pair density calculated in the approximate 
analytic fashion. It was not pursued because the analytic method 
proved more fruitful, and convergence problems arose in the numeri- 
cal method at large optical depths. In both cases, though, no annihil- 
ation feature was determined to be expected from a stable, relativistic 
plasma in pair balance, agreeing with the conclusion of Svensson (1984). 
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Figure 13. Numerical simulation of nonmagnetic relativistic 

thermal plasma. Production and emergent spectra of a 
uniform nonmagnetic relativistic plasma are shown for 
temperatures and optical depths indicated. Calculated 
equilibrium positron densities are also shown. Label "1" 
refers to the photon production spectrum from bremsstrah- 
lung and annihilation processes, label "2" refers to the 
steady photon spectrum after suffering effects of absorption, 
and label "3" refers to the emergent spectrum after 
Comptonization. 
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V. SUMMARY AND CONCLUSIONS 

This preliminary investigation into the physics of nonmag- 
netic, relativistic thermal plasmas was motivated by astrophysical 
observations showing that the existence of plasmas at relativistic 
temperatures may not be an uncommon feature of the universe. Most 
models suggest that relativistic plasmas originate from the dissipa- 
tion of kinetic or elastic energy in the intense gravitational field of 
a compact massive object. Other speculations, such as bulk matter- 
antimatter annihilation, must however not be excluded. 

In order to frame quantitative analyses of models relating to 
such observations, the elementary physics must be properly per- 
formed. It was the object of this study to help contribute to this 
goal. The principal results of this work are as follows: 

(1) A r elativistically correct description of the kinematics 
of particle and photon gases was performed, giving exact expres- 
sions for reaction rates and opacities of an isotropic gas consisting 
of photons and/or particles, and luminosities and production spectra 
from an isotropic particle plasma. 

(2) Specializing to an isotropic Maxwell- Boltzmann particle 
gas consisting generally of different mas s particles led to expres- 
sions for the reaction rate and luminosity in terms of a single inte- 
gral over the invariant cross section, and to the production spectra 
(yielding either photons or material particles) in terms of a double 
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integral over the center -of-momentum cross section differential 
in the energy of the produced particles or photons. 

(3) These results were used to analytically calculate the reac- 
tion rate, luminosity, and production spectra of annihilation photons 
from a Maxwell- Boltzmann electron- positron plasma, and to numeri- 
cally calculate thermal electron-electron bremsstrahlung . Thermal 
electron- positron bremsstrahlung was calculated by interpolating 
between the modified Bethe -Heitler formula at nonrelativistic ener- 
gies and the lowest order electron-electron bremsstrahlung cross 
section at relativistic energies. 

(4) A treatment of the Compton power and moments of the 
Compton scattered spectrum was used to provide an approximate 
description of Comptonization in a relativistic plasma of photons and 
particles. 

(5) Single particle reaction rates and energy loss timescales 
were examined in a relativistic ally correct manner. An expression 
for the luminosity (energy loss or energy exchange rate per unit 
volume) from the interaction of two different temperatures , different 
mass Maxwell- Boltzmann gases of particles was derived. This has 
potential application to two-temperature accretion disk models. 

(6) Calculation of energy loss timescales led to the conclu- 
sion that an optically thin relativistic electron-proton plasma cannot 
obtain the Maxwell- Boltzmann form at temperatures greater than 
about 1 MeV, agreeing with the results of other workers. 



(7) A compilation of pion production data from proton-proton 
collisions was made from threshold to ISR energies. The reaction 
rates and luminosities for different modes of pion production were 
calculated for a Maxwell- Boltzmann gas of protons at temperatures 
between 1 MeV and 1 GeV. It was concluded that above temperatures 
greater than about 60 MeV protons in the gas cannot obtain the 
Maxwell- Boltzmann form because of the energy loss due to pion 
production. 

(8) A consideration of dominant processes was used to model 
a steady, uniform, nonmagnetic, relativistic thermal plasma. An 
approximate treatment of Comptonization was developed at subrela- 
tivistic temperatures 0^0.4 and extreme relativistic tempera- 
tures 0^3, where 0 is the temperature measured in units of the 
rest mass energy of the electron. Asymptotic expressions for rate 
processes were calculated, permitting quantitative analysis of the 
system. 

(9) The equilibrium pair density was calculated as a function 
of temperature, proton optical depth, and proton number density. 
Either zero, one, or two roots were found for each set of values. 
Above a critical temperature dependent on the proton optical depth 
and proton number density, only a system in full thermodynamic 
equilibrium could exist under conditions of steady state. 
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(10) No annihilation feature should be observed from a non- 
magnetic relativistic thermal plasma in pair equilibrium, in agree- 
ment with other studies. 

Further work must concern the properties of relativistic 
thermal plasmas in the presence of magnetic fields. The first steps 
in this direction have been taken at extreme relativistic tempera- 
tures (Araki and Suguru 1983; Kusunose and Takahara 1983). Having 
characterized the static properties of these systems, a dynamic 
treatment will be possible, permitting a better understanding of such 
interesting systems as gamma ray bursts, active galaxies, extra- 
galactic jets, and quasi- stellar objects. 



APPENDIX: NUMERICAL CALCULATION OF THERMAL 
BREMSSTRAHLUNG 

Calculation of the Gaunt factors associated with thermal e-p, 
e-e, and e + -e brems strahlung was performed numerically on a 
VAX/VMS. For e-p brems strahlung, only a single integral had to 
be evaluated (equation (III-14)). The integration was performed 
using Gaus s-Laguerre quadrature after making the variable change 
y = (y - 1 -k). A 20-point routine proved sufficient. 

The three-fold integral for the calculation of thermal e-e 
brems strahlung, equation (III- 16), was more difficult due to the peak- 
ing of the angle -dependent differential cross section in the forward 
and backward directions. The small angle regime was magnified 
through the transformation cos 0 = 1- exp(-z), extending the range 
of the integral from 0 to 00 Also, double precision was required 
in the angular integration because of the presence of large, nearly 
cancelling factors. 

For photon energies that were not small, that is, for 

k ^ © /10, the transformation x = ? /(k ; + 1) served to broaden the 

e ee 

range of the k / integration. In both cases, the 20-point Gauss- 
Laguerre quadrature formula was employed for the integration. A 
Simpson algorithm, written to evaluate an integral to any desired 
degree of accuracy, was used to perform the integration over y r • 
Because of the slow variation of the integrand in the variable y r » 

8 or 16 iterations were usually enough to calculate the integral over 
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an interval in to 1% accuracy, and typically 5 to 10 intervals of 

y were needed to calculate the Gaunt factor to the same accuracy. 
This method exhibited which relative Lorentz factors are chiefly 
responsible for bremsstrahlung emission at a given temperature and 
photon energy. 

For extremely low photon energies, the Gauss-Laguerre 
technique was not sufficiently accurate in the integration over k* 
due to the effects of the exponential in equation (III- 18). This caused 
the value of the integral to be due largely from contributions in the 
immediate neighborhood of k' = k. In this case, the integral over 
k* was evaluated with Simpson's approximation to either side of 
k' = k, after determining in the program the range of the exponential 
factor for each iteration. 

•f- — 

The Gaunt factors for thermal e -e bremsstrahlung were 

evaluated as for the e-e case, with a cross section appropriate to 

4 * - 

the problem. At relativistic energies, the spectra of e -e and 
e-e bremsstrahlung are the same except at the hard ends of the 
spectra (Baier, Fadin, and Khoze 1968), which are not supposed to 
contribute appreciably after averaging over a thermal distribution 
(Svensson 1982). Thus the Gaunt factor for e"^-e bremsstrahlung 
is a factor of 2 grester than e-e bremsstrahlung at relativistic tem- 
peratures as the electron and positon are distinct. 

At nonrelativistic (NR) energies, the bremsstrahlung emis- 
sion in e + -e scattering is larger than the e-e emission since the 
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e -e system possesses a dipole moment. Joseph and Rohrlich 
(1958) argue that the CM bremsstrahlung spectrum from e + -e scat- 
tering in a collision characterized by the relative velocity 8^ is 
equivalent to the emission from a charged particle of reduced mass 
m e /2 scattering with an infinitely heavy proton, again at relative 
velocity 8^ . This is simply the NR B-H formula with m^ -> m^/2: 


f + J k ’ Y r ) 

e e 


a r 


dCJ B-H (k ’ ^ r> 
dk 


m -> m / 2 
e e 


y - 1 « 1 

' r 


16 

3P 2 


ton 


8 + 8 ' 


(Al) 


/ / / 2 
The value of |3 is defined through y = y - 2k. Hence 8 = 

O 2 , 

p - 4k in the NR limit (whereas y = y^ - k for e-p bremsstrah- 
lung). 

On this basis, the low energy ( y^ < 2) e -e bremsstrahlung 

was taken to be the B-H formula, with reduced mass m /2. In the 

e 

trans relativistic regime, as shown in Figure A, the B-H spectrum 

(solid line) and the e-e bremsstrahlung spectrum (dotted line), given 

there as a function of UU = m e ^» are found to be agreeably similar. 

For 2 < y , ^ 4, the e + -e bremsstrahlung spectrum is taken to be 

l/2(y - 2)f + 1/2(4- y ) f TT , as shown explicitly for the case y 

= 3 by the dot-dashed line. At the two extremes, the appropriate 

spectra are used as discussed above. Also, the exact value § = 

e-e 
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1 / 2 x 

(y - l)/[2(y +1)] for the maximum CM photon energy is used in 
r r 

the interpolation region. At NR energies, l/2(y - 1) -» ? , as 

r e - e 

required. 

At this point, I wish to correct a mistake made by Svensson 

(1982). He claims that the NR e"*"-e bremsstrahlung luminosity 

3/2 

evaluated through the NR cross section (Al) is 2 greater than the 

_ 1/2 

NR e-p bremsstrahlung luminosity. Actually, this ratio is 2 
First note that in the limit y . - 1 « 1, the integral 


e (Y r )a(v r ) 


1 f e ‘ e .„ , d °B-H (k; V . 

2 I dkk dk v - 1 « 1 

a r r 

e 0 


, / f P r /4 , ,/B + |3' 

T I dk ~2 ^(fF 7 

P r 

0 r 


can be performed exactly and is independent of |3 . Its value is 1 

(in units of the electron rest mass energy) for e-p scattering and 
1/2 for e + -e scattering. Therefore, the ratio of luminosities, 
other than this factor of 2, simply reflects the relative number of 
scatterings per unit volume per unit time (the reaction rate) in the 
following two systems: (i) NR electrons scattering with NR protons, 

and (ii) NR electrons scattering with NR positrons, all scatterings 
defined by the same energy-independent cross section. The ratio 


of the reaction rates of (ii) to (i) is just 2 , as can be seen either 
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through equations (11-16) or (11-28), or simple considerations of the 
NR reaction rates. Since each e + -e scattering yields, on average, 
one-half the energy as an e-p scattering, but occurs at a rate greater 
by a factor of 2^^ , the conclusion follows. 
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Figure A. Single particle bremsstrahlung spectra. Frequency- 

dependent bremsstrahlung spectra as a function of the 

relative Lorentz factor v are shown for different values 

r 

of y in the interpolation regime for the approximate 
calculation of thermal e + -e bremsstrahlung. The solid 
line refers to the modified Bethe-Heitler formula, employ- 
ing the reduced mass m /2. The dotted line refers to the 

e 

exact, lowest order e -e bremsstrahlung spectrum, and 
the interpolation spectrum, shown only for = 3, is 
given by the dot- da shed line. The formula for the inter- 
polation spectrum is given in the text. 
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